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A  theoretical  research  program  directed  toward  the  study  of  the  energe¬ 
tics  and  LWIR  radiative  properties  of  selected  uranium/ oxygen  band  systems  has 
been  undertaken.  Included  in  this  research  program  was  the  investigation  of 
the  strongest  electronic  and  vibrational  bands  in  the  LWIR  region  for  the 
species  UO,  UO'*',  U02»  U02^.  and  U02*'*’*  The  program  for  accomplishing  this 
research  effort  was  formulated  into  three  separate  tasks:  a)  adaption  of  our 
electronic  structure  codes  to  the  DNA  CYBER  176  System,  b)  calculation  of  per¬ 
tinent  electronic  wavefunctions  and  energies,  as  a  function  of  internuclear 
separation  and  within  a  relativistic  framework,  for  selected  species  of  the 
uranium/oxygen  system  which  may  be  important  in  the  LWIR  region,  and  c)  cal¬ 
culation  of  electronic  transition  moments  and  transition  probabilities  between 
specific  vibrational  levels  of  the  electronic  states  corresponding  to  the 
strongest  radiating  band  systems  belonging  to  the  uranium/ oxygen  system  and 
prediction  of  IR  and  possible  optical  oscillator  strengths. 

Our  calculations  indicate  that  the  species  UO'*'  will  be  efficiently  solar 
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~  33000  cm“ ‘  IS  the  first  strongly  optically  connected  state.  The  wavelength 

for  excitation  of  this  and  higher  coupled  electronically  excited  states  is 

A  <  300  nra,  which  places  them  beyond  the  main  region  of  the  solar  flux.  The 

neutral  UO 2  species  has  the  first  strongly  coupled  optical  transition  at 

A  ~  2.60  nra,  again  beyond  the  region  for  efficient  solar  pumping.  The 


pathways  for 'convers ion  of  solar  photons  to  the  LWIR  region  are  still 
uncertain  but  strong  LWIR  radiation  (fgj  ~  1.2  x  10““*)  is  predicted  for  the 
vibrational  transitions  of  the  ground  state  of  UO-'*’.  Further  studies  of 
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SECTION  1 


INTRODUCTION 

The  release  of  certain  chemical  species  into  the  upper  atmosphere  results 
in  luminous  clouds  that  display  the  resonance  electronic-vibration-rotation 
spectra  of  the  chemically  reacting  species.  Such  spectra  are  seen  in  rocket 
releases  of  chemicals  for  upper  atmospheric  studies,  upon  re-entry  into  the 
atmosphere  of  artificial  satellites  and  missiles,  and  as  a  result  of  energy 
deposition  in  the  atmosphere  caused  by  nuclear  weapons  effects.  Of  particular 
interest  in  this  connection  is  the  observed  spectra  of  certain  metallic 
oxides.  Prom  band  intensity  distributions  of  the  spectra,  and  knowledge  of 
the  f-values  for  electronic  and  vibrational  transitions,  the  local  conditions 
of  the  atmosphere  can  be  determined  (Reference  1).  Such  data  are  fundamental 
for  the  analysis  of  detection  and  discrimination  problems. 

Present  theoretical  efforts,  which  are  directed  toward  a  more  complete 
and  realistic  analysis  of  the  transport  equations  governing  atmospheric 
relaxation  and  the  propagation  of  artificial  disturbances,  require  detailed 
information  of  thermal  opacities  and  LWIR  absorption  in  region  of  temperature 
and  pressure  where  both  atomic  and  molecular  effects  are  important  (References 
2  and  3).  Although  various  experimental  techniques  have  been  employed  for 
both  atomic  and  molecular  systems,  theoretical  studies  have  been  largely 
confined  to  an  analysis  of  the  properties  (bound-bound,  bound-free  and  free- 
free)  of  atomic  systems  (References  4  and  5).  This  has  been  due  in  large  part 
to  the  unavailability  of  reliable  wavefunctions  for  diatomic  molecular 
systems,  and  particularly  for  excited  states  or  states  of  open-shell 
structure.  Only  recently  (References  6-8)  have  reliable  procedures  been 
prescribed  for  such  systems  which  have  resulted  in  the  development  of 
practical  computational  programs. 

The  application  of  these  computational  methods  to  studies  of  the  elec¬ 
tronic  structure  and  radiation  characteristics  of  metal  oxides  has  been 
reported  for  several  of  the  lighter  systems  (References  9-11).  A  preliminary 
study  of  the  uranium/oxygen  system  has  been  reported  by  Michels  (Reference 
12)  which  identified  a  large  number  of  low-lying  molecular  states  for  both  the 


UO  and  UO'*’  systems.  Of  particular  interest  was  the  discovery  of  two  struc¬ 
tures  for  UO'*'  that  resulted  from  two  different  spin-couplings  of  the  uranium 
valence  electrons.  These  results  suggested  strong  LWIR  radiation  in  UO'*’  aris¬ 
ing  from  different  electronic  state  transitions. 

An  inherent  uncertainty  in  these  preliminary  calculations  was  present, 
owing  to  the  neglect  of  relativistic  effects  that  were  much  too  difficult  to 
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include  in  molecular  calculations  at  that  period  of  time.  The  7s  valence 
electron  of  uranium,  and  its  corresponding  o-bonding  molecular  orbital,  are 
highly  relativistic  in  nature  which  results  in  a  contracted  charge  density 
relative  to  that  which  would  occur  in  lighter  molecular  systems.  The  effect 
of  this  contraction  on  the  relative  positions  of  the  low-lying  electronic 
states  of  the  uranium/oxygen  system  can  now  be  calculated  with  some  degree  of 
confidence  using  newly  developed  relativistic  computer  codes. 

Because  of  inherent  difficulties  in  the  experimental  determination  of  the 
spectroscopy,  transition  probabilities  and  LWIR  radiation  for  metal  oxide 
systems  and  in  light  of  the  aforementioned  recent  progress  in  the  calculation 
of  relativistic  electronic  wavefunct ions ,  especially  for  diatomic  systems,  a 
technical  program  for  calculating  these  properties  was  undertaken  for  the 
Defense  Nuclear  Agency  under  Contract  DNA001-82-C-0015.  The  emphasis  in  this 
work  was  on  the  ions  of  uranium  and  uranium  oxide  (U^,  UO^)  since  these  have 
been  determined  to  be  important  radiators  in  the  LWIR  region.  These  studies 
indicated  that  80^^  and  the  doubly  ionized  species,  and  UO^^,  should 

also  be  considered  because  of  their  role  in  charge  neutralization  processes 
and  their  potential  as  early-time  radiators. 

A  careful  re-evaluation  of  the  ground  and  low-lying  excited  states  of 
UO^  was  carried  out  under  the  current  research  program.  Earlier  studies  have 
suggested  the  possibility  of  strong  absorption  for  this  species  in  the  solar 
pumped  wavelength  region.  More  definitive  studies  which  rule  against  this 
possibility  are  presented  in  this  report. 

The  general  composition  of  this  report  is  as  follows.  In  Section  II,  we 
present  a  description  of  the  mathematical  methods  which  were  employed  in  this 
research.  Included  in  Section  II  are  sub-sections  which  deal  with  the 
construction  of  electronic  wavefunct ions ,  the  calculations  of  expectation 
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properties,  the  evaluation  of  molecular  transition  probabilities,  and  the 
calculation  of  electronic  wavefunctions  using  both  the  ab  initio  and  density 
functional  methods.  This  is  followed  by  Section  III  which  describes  the 
inclusion  of  relativistic  effects  into  the  density  functional  (X^)  method. 
The  calculated  results  and  pertinent  discussion  are  presented  in  Section  IV. 
Recommendations  are  presented  in  Section  V. 


SECTION  2 


METHOD  OF  APPROACH  -  NONRELATIVISTIC  METHODS 

2.1  QUANTUM  MECHANICAL  CALCULATIONS 

Central  to  these  theoretical  studies  are  the  actual  quantum-mechanical 
calculations  which  must  be  carried  out  for  the  atomic  and  molecular  species. 
For  added  clarity,  various  aspects  of  these  calculations  are  discussed  in 
individual  subsections. 

2.1.1  Levels  of  Approximation 

Much  evidence  on  diatomic  and  polyatomic  systems  indicates  the  inadequacy 
of  a  minimum  Slater-type-orbital  (STO)  basis  for  constructing  quantitatively 
correct  molecular  wavefunct ions  (References  13  and  14).  This  means  inner- 
shell  and  valence-shell  STO's  of  quantum  numbers  appropriate  to  the  atoms  (Is, 
2s,  2p,  for  C,  N,  0;  etc.).  The  main  deficiency  of  the  minimum  basis  set  is 
in  its  inability  to  properly  describe  polarization  and  the  change  of  orbital 
shape  for  systems  which  exhibit  large  charge  transfer  effects.  Values  of  the 
screening  parameters  C  for  each  orbital  can  either  be  set  from  atomic  studies 
or  optimized  in  the  molecule;  the  latter  approach  is  indicated  for  studies  of 
higher  precision.  When  high  chemical  accuracy  is  required,  as  for  the 
detailed  studies  of  the  ground  or  a  particular  excited  state  of  a  system,  a 
more  extended  basis  must  be  used.  Double-zeta  plus  polarization  functions  or 
optimized  MO's  are  required  for  reliable  calculated  results  of  chemical 
accuracy . 

The  chosen  basis  sets  give  good  results  only  when  used  in  a  maximally 
flexible  manner.  This  implies  the  construction  of  Cl  wavefunctions  with  all 
kinds  of  possible  orbital  occupancies,  so  that  the  correlation  of  electrons 
into  overall  states  can  adjust  to  an  optimum  form  at  each  geometrical 
conformation  and  for  each  state.  Except  when  well-defined  pairings  exist,  as 
for  closed  shell  and  exchange  dominated  systems,  a  single-configuration  study 
(even  of  Hartree-Fock  quality)  will  be  inadequate. 

2.1.2  Spin  and  Symmetry 

Proper  electronic  states  for  systems  composed  of  light  atoms  should 
possess  definite  eigenvalues  of  the  spin  operator  S^  as  well  as  an  appropriate 
geometrical  symmetry.  The  geometrical  symmetry  can  be  controlled  by  the 
assignment  of  orbitals  to  each  configuration,  but  the  spin  state  must  be 
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obtained  by  a  constructive  or  projective  technique.  Formulas  have  been 
developed  (Reference  15)  for  projected  construction  of  spin  states  from 
orthogonal  orbitals,  and  programs  implementing  these  formulas  have  been  in 
routine  use  at  UTRC  for  several  years. 

One  of  the  least  widely  appreciated  aspects  of  the  spin-projection  prob¬ 
lem  is  that  the  same  set  of  occupied  spatial  orbitals  can  sometimes  be  coupled 
to  give  more  than  one  overall  state  of  given  S  quantum  number.  It  is  neces¬ 
sary  to  include  in  calculations  all  such  spin  couplings,  as  the  optimum  coup¬ 
ling  will  continuously  change  with  changes  in  the  molecular  conformation. 

This  is  especially  important  in  describing  degenerate  or  near-degenerate 
excited  electronic  states. 

2.1.3  Method  of  Ab  Initio  Calculation 

A  spin-free,  nonrelat ivist ic ,  electrostatic  Hamiltonian  is  employed  in 
the  Born-Oppenheimer  approximation.  In  systems  containing  atoms  as  heavy  as 
Kr,  this  approximation  is  quite  good  for  low-lying  molecular  states.  For  a 
diatomic  molecule  containing  n  electrons,  the  approximation  leads  to  an  elec¬ 
trostatic  Hamiltonian  depending  parametrically  on  the  internuclear  separation, 
R; 
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where  and  Zg  are  the  charges  of  nuclei  A  and  B,  and  r^  is  the  separation 
of  electron  i  and  nucleus  A.  ^  is  in  atomic  units  (energy  in  hartrees,  length 
in  bohrs) . 

Electronic  wavefunct ions  tiR)  are  made  to  be  optimum  approximations  to 
solutions,  for  a  given  R,  of  the  Schrodinger  equation 


<^(R)'KR)=  E  (R)'I'(R) 

by  invoking  the  variational  principle 

/'i'*(R)«^{R)'KR)(^T 
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The  integrations  in  Equation  (3)  are  over  all  electronic  coordinates  and 
the  stationary  values  of  W(R)  are  approximations  to  the  energies  of  states 
described  by  the  corresponding  \KR)<  States  of  a  particular  symmetry  are 
studied  by  restricting  the  electronic  wavefunction  to  be  a  projection  of  the 
appropriate  angular  momentum  and  spin  operators.  Excited  electronic  states 
corresponding  to  a  particular  symmetry  are  handled  by  construction  of  config¬ 
uration-interaction  wavefunctions  of  appropriate  size  and  form. 

The  specific  form  for  i(/(R)  may  be  written 


where  each  is  referred  to  as  a  configuration,  and  has  the  general  struc- 


where  each  (Ji^  is  a  spatial  orbital,  /7  is  the  antisymmetrizing  operator, 

is  the  spin-projection  operator  for  spin  quantum  number  S,  and  is  a  product 

of  a  and  S  one-electron  spin  functions  of  magnetic  quantum  number  M.  No 

requirement  is  imposed  as  to  the  double  occupancy  of  the  spatial  orbital,  so 

Equations  (4)  and  (5)  can  describe  a  completely  general  wavefunction. 

In  Hartree-Fock  calculations  i|<R)  is  restricted  to  a  single  which  is 

assumed  to  consist  as  nearly  as  possible  of  doubly-occupied  orbitals.  The 

orbitals  .  are  then  selected  to  be  the  linear  combinations  of  basis  orbitals 
W. 

best  satisfying  Equation  (3).  Writing 

4>fj.\  =  T  X  V  (6) 

the  a^£  are  determined  by  solving  the  matrix  Hartree-Fock  equations 

I  (each  X) 

where  is  the  orbital  energy  of 


(7) 


The  Fock  operator  has  been  thoroughly  discussed  in  the  literature 
(Reference  16)  and  depends  upon  one-  and  two-electron  molecular  integrals  and 
upon  the  This  makes  Equation  (7)  nonlinear  and  it  is  therefore  solved 

iteratively.  UTRC  has  developed  programs  for  solving  Equation  (7)  for  both 
closed  and  open-shell  systems,  using  basis  sets  consisting  of  Slater-type 
atomic  orbitals.  Examples  of  their  use  are  in  the  literature  (Reference 
7). 

In  configuration  interaction  calculations,  the  summation  of  Equation  (4) 
has  more  than  one  term,  and  the  c^  are  determined  by  imposing  Equation  (3)  to 
obtain  the  secular  equation 

Z  {Huv“  wSuv)c^  =  0  (eochu.)  (8) 
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Equation  (8)  is  solved  by  matrix  diagonal izat ion  using  either  a  modified 
Givens  method  (Reference  17)  or  a  method  due  to  Shavitt  (Reference  18). 

The  matrix  elements  and  S^^may  be  reduced  by  appropriate  operator 

algebra  to  the  forms 
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where  P  is  a  permutation  and  e  its  parity.  The  sum  is  over  all  permutations. 

P 

<^|  gP|  is  a  "Sanibel  coefficient"  and  the  remaining  factors  are  spatial 

integrals  which  can  be  factored  into  one-  and  two-electron  integrals.  If  the 

d)  .  are  orthonormal,  Equations  (10)  and  (11)  become  more  tractable  and  the 

m 

H  and  S  may  be  evaluated  by  explicit  methods  given  in  the  literature 


(Reference  15).  Computer  programs  have  been  developed  for  carrying  out  this 
procedure,  and  they  have  been  used  for  problems  containing  up  to  106  total 
electrons,  10  unpaired  electrons,  and  several  thousand  configurations. 

The  Cl  studies  described  above  can  be  carried  out  for  any  orthonormal  set 
of  for  which  the  molecular  integrals  can  be  calculated.  Programs  devel¬ 
oped  by  UTRC  make  specific  provision  for  the  choice  of  the  as  Slater-type 
atomic  orbitals,  as  symmetry  molecular  orbitals,  as  Hartree-Fock  orbitals,  or 
as  more  arbitrary  combinations  of  atomic  orbitals. 

2.1.4  Molecular  Integrals 

The  one-  and  two-electron  integrals  needed  for  the  above  described  method 
of  calculation  are  evaluated  for  STO's  by  methods  developed  by  the  present 
investigators  (Reference  19).  All  needed  computer  programs  have  been  devel¬ 
oped  and  fully  tested  at  UTRC. 

2.1.5  Configuration  Selection 

Using  a  minimum  basis  plus  polarization  set  of  one-electron  functions,  a 
typical  system  can  have  of  the  order  of  10**  configurations  in  full  Cl  (that 
resulting  from  all  possible  orbital  occupancies).  It  is  therefore  essential 
to  identify  and  use  the  configurations  describing  the  significant  part  of  the 
wavefunction.  There  are  several  ways  to  accomplish  this  objective.  First, 
one  may  screen  atomic-orbital  occupancies  to  eliminate  configurations  with 
excessive  numbers  of  anti-bonding  orbitals.  A  third  possibility  is  to  carry 
out  an  initial  screening  of  configurations,  rejecting  those  whose  diagonal 
energies  and  interaction  matrix  elements  do  not  satisfy  significance  criteria. 
Programs  to  sort  configurations  on  all  Che  above  criteria  are  available  at 
UTRC. 

Other,  potentially  more  elegant  methods  of  configuration  choice  involve 
formal  approaches  based  on  natural-orbital  (Reference  20)  or  multiconfigura¬ 
tion  SCF  (Reference  21)  concepts.  To  implement  the  natural-orbital  approach, 
an  initial  limited-CI  wavefunction  is  transformed  to  natural-orbital  form,  and 
the  resulting  natural  orbitals  are  used  to  form  a  new  Cl.  The  hoped-for 
result  is  a  concentration  of  the  bulk  of  the  Cl  wavefunction  into  a  smaller 
number  of  significant  terms.  The  multiconfiguration  SCF  approach  is  more 
cumbersome,  but  in  principle  more  effective.  It  yields  the  optimum  orbital 
choice  for  a  preselected  set  of  configurations.  This  approach  works  well  when 
a  small  number  of  dominant  configurations  can  be  readily  identified. 


It  should  be  emphasized  that  the  problem  of  configuration  choice  is  not 
trivial,  and  represents  an  area  of  detailed  study  in  this  research.  The  exis¬ 
tence  of  this  problem  causes  integral  evaluation  to  be  far  from  a  unique 
limiting  factor  in  the  work. 

2.1.6  Density  Function  Approach  -  Xa  Model 

The  Xa  model  (Reference  22)  for  the  electronic  structure  of  atoms,  mole¬ 
cules,  clusters  and  solids  is  a  local  potential  model  obtained  by  making  a 
simple  approximation  to  the  exchange  -  correlation  energy.  If  we  assume  a 
nonrelativistic  Hamiltonian  with  only  electrostatic  interactions,  it  can  be 
shown  that  the  total  energy  E  of  a  system  can  be  written  exactly  (Reference 
23)  (in  atomic  units)  as 


E=  z  n,  <  Ui  --5- 7,  +Z  -7^  U.  >  +  4“  5!  -j 
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This  expression  is  exact  provided  the  u£  are  natural  orbitals  and  n^  are 
their  occupation  numbers  (i.e.,  eigenfunctions  and  eigenvalues  of  the  first 
order  density  matrix).  The  first  term  in  Equation  (12)  represents  the  kinetic 
and  electron-nuclear  energies.  The  second  term  is  the  nuclear  repulsion 
energy.  The  sums  (u,v)  are  over  all  the  nuclear  charges  in  the  system.  The 
third  term  is  the  electron-electron  repulsion  term,  which  represents  the  clas¬ 
sical  electrostatic  energy  of  the  charge  density  p  interacting  with  itself. 


where 


^(1)  =Z  n.u*(l)  Uj(i) 


The  last  term  E^^  represents  the  exchange  correlation  energy  and  can  be 
expressed  formally  as 


=  d'T/ -Sssiiii  ar,  , 


where  (1,  2)  represents  the  exchange-correlation  hole  around  an  electron 

at  position  1.  In  the  exact  expression,  P  is  dependent  on  the 
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secoad-order  density  matrix.  In  the  Hartree-Fock  approximation  is  the 
exchange  energy,  represents  the  Fermi  hole  due  to  the  exclusion  principle 

and  depends  only  on  the  first-order  density  matrix.  In  the  Xa  method,  we  make 
a  simpler  assumption  about  if  wfi  assume  that  the  exchange-correlation 

hole  is  centered  on  the  electron  and  is  spherically  symmetric,  it  can  be  shown 
that  the  exchange-correlation  potential 


U  =  f  • 
xc  J 


is  inversely  proportional  to  the  range  of  the  hole,  rg ,  where  rg  is 
defined  by 

47r 

~ir~  r}  /){0  =  l 


Therefore,  in  the  Xa  model,  the  potential  is  proportional  to  p  (r).  We 


define  a  scaling  parameter  a  such  that 


The  expression  in  Equation  (37)  is  defined  so  that  a  =  2/3  for  the  case  of  a 
free  electron  gas  in  the  Hartree-Fock  model  (Reference  24)  and  a  =  1  for  the 
potential  originally  suggested  by  Slater  (Reference  25).  A  convenient  way  to 
choose  this  parameter  for  molecular  and  solid  state  applications  is  to  opti¬ 
mize  the  solutions  to  the  Xo  equations  in  the  atomic  limit.  Schwarz  (Refer¬ 
ence  26)  has  done  this  for  atoms  from  z=  1  to  z=41  and  found  values  between 
2/3  and  1. 

In  the  "spin  polarized"  version  of  the  Xa  theory,  it  is  assumed  (as  in 
the  spin-unrestricted  Hartree-Fock  model)  that  electrons  interact  only  with  a 
potential  determined  by  the  charge  density  of  the  same  spin.  In  this  case  the 
contribution  to  the  total  energy  is  summed  over  the  two  spins,  s  =  ±  1/2. 
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where  the  potential  i.  spin-dependent 
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and  pg  is  the  charge  density  corresponding  to  electrons  of  spin  s.  The 
spin  polarized  Xa  model  is  useful  for  describing  atoms  and  molecules  with 
open-shell  configurations  and  crystals  which  are  ferromagnetic  or  anti- 
ferromagnet ic . 

Once  one  has  made  the  Xa  approximation  to  the  total  energy  functional  E 
in  Equation  (12),  then  the  rest  of  the  theory  follows  from  the  application  of 
the  variational  principle.  The  orbitals  uj_  are  determined  by  demanding  that 
E  be  stationary  with  respect  to  variations  in  u£.  This  leads  to  the  set  of 
one-electron  X a  equations 


where  is  the  one-electron  eigenvalue  associated  with  u^.  Since  p(r)  is 
defined  in  terras  of  the  orbitals  u^,  Equation  (20)  must  be  solved  itera¬ 
tively,  until  self-consistency  is  achieved.  Empirically,  if  one  takes  as  an 
initial  guess  that  p  is  approximately  a  sum  of  superimposed  atomic  charge 
densities,  then  the  convergence  of  this  procedure  is  fairly  rapid.  The  factor 
of  2/3  multiplying  the  potential  is  a  result  of  the  linear  dependence  of  E^^ 
on  p.  This  also  has  a  consequence  that  the  Xa  eigenvalues  do  not  satisfy 
Koopraan's  theorem,  i.e.,  they  cannot  be  interpreted  as  ionization  energies. 
However,  it  can  be  shown  that  the  Cj  are  partial  derivatives  of  the  total 
expression  of  Equation  (12)  with  respect  to  the  occupation  number. 


If  E  were  a  linear  function  of  n^,  then  Koopmans*  theorem  would  hold. 

However,  because  of  the  dominant  Coulomb  term,  E  is  better  approximated  by  a 
quadratic  function  in  n£ .  This  leads  to  the  "transition  state"  approxima¬ 
tion  which  allows  one  to  equate  the  difference  in  total  energy  between  the 
state  (n.,  n.)  and  (n.-  1,  n.+l)  to  the  difference  in  the  one-electron 

i’  j  1  j 

energies  e.  -  e.  calculated  in  the  state  (n.  -  1/2,  n.  +  1/2).  The  error  in 
J  i  1-  J 

this  approximation  is  proportional  to  third-order  derivatives  of  E  with 

respect  to  n^  and  n j ,  which  are  usually  small  (Reference  27).  The  main 

advantage  of  using  the  transition  state  rather  than  directly  comparing  the 

total  energy  values  is  computational  convenience,  especially  if  the  total 

energies  are  large  numbers  and  the  difference  is  small. 

The  relationship  of  Equation  (21)  also  implies  the  existence  of  a  "Fermi 

level"  for  the  ground  state.  This  can  be  seen  by  varying  E  with  respect  to 

E 

n£  under  the  condition  that  the  sum  j  is  a  constant,  i.e.. 


(22) 


implies  3E/ 3nj^  =  X,  where  X  is  a  Lagrangian  multiplier.  This  implies 
the  total  energy  is  stationary  when  all  the  one-electron  energies  are 
However,  the  occupation  numbers  are  also  subject  to  the  restriction  0 
This  leads  to  the  following  conditions  on  the  ground  state  occupation 
numbers ; 

«:  <X)  n. =1 

*  I 

€j  >X  )  n.  =0 

€i=X  )•  0<nj<i 


that 
equal . 

<  ni  ^  1 . 
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In  other  words,  the  ground  state  eigenvalues  obey  Fermi  statistics  with  X 
representing  the  Fermi  energy.  It  should  be  noted  that,  in  contrast  to  the 
Hartree-Fock  theory,  where  all  the  n^  are  either  0  or  1,  the  Xa  model  pre¬ 
dicts,  in  some  cases,  fractional  occupation  numbers  at  the  Fermi  level.  In 
particular,  this  will  occur  in  a  system  (such  as  transition  metal  or  actinide 
atom)  which  has  more  than  one  open  shell. 
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The  Xa  model  differs  in  other  significant  ways  from  the  Hartree-Fock 
method.  In  fact,  the  simplification  introduced  in  approximating  the  total 
energy  expression  introduces  several  distinct  advantages  over  Hartree- 
Fock: 

1.  The  primary  advantage  is  purely  computational.  The  one-electron 
potential  in  Equation  (20)  is  orbital-independent  and  local,  i.e.,  it  is  the 
same  for  all  electrons  (except  in  the  spin-polarized  Xa  theory)  and  is  a 
multiplicative  operator.  On  the  other  hand,  the  Hartree-Fock  potential  is 
nonlocal,  or  equivalently,  there  is  a  different  local  potential  for  each 
orbital.  This  involves  a  great  deal  more  computational  effort,  especially  for 
systems  described  by  a  large  number  of  orbitals.  It  has  been  shown  (Reference 
28)  that  the  Xa  orbitals  for  the  first  and  second  row  atoms  are  at  least  as 
accurate  as  a  double-zeta  basis  set,  and  are  probably  better  for  larger  atoms 
which  involve  electrons  with  A  >  2. 

2.  The  orbital-independent  Xa  poten:ial  leads  to  a  better  one-electron 
description  of  electronic  excitations  of  a  system.  Both  the  unoccupied 

(n.  =  0)  and  occupied  (n.  =  1)  eigenfunctions  are  under  the  influence  of  the 
same  potential  resulting  from  the  other  N-1  electrons.  The  Hartree-Fock  vir¬ 
tual  orbitals  see  a  potential  characteristic  of  the  N  occupied  orbitals,  and 
therefore  are  not  as  suitable  for  describing  the  excited  states.  Actually, 
although  the  ground  state  virtual  eigenvalues  are  usually  a  good  description 
of  the  one-electron  excitations,  the  virtual  spectrum  of  the  transition  state 
potential  where  one-half  an  electron  has  been  removed  from  the  system  gives  a 
much  better  first-order  picture  of  these  levels  (Reference  29). 

3.  As  has  been  shown  by  Slater  (Reference  30),  the  Xa  model  rigorously 
satisfies  both  the  virtual  and  Hel Iman-Feynman  theorems,  independent  of  the 
value  of  the  parameter  a.  This  is  convenient  for  calculating  the  force  on  a 
nucleus  directly  in  terms  of  a  three-dimensional  integral,  rather  than  the 
six-dimensional  integrals  in  the  expression  for  the  total  energy  of  Equation 
(12)  . 

2.1.7  Computational  Aspects  of  the  Xa  Method 

In  application  of  the  Xa  model  to  finite  molecular  systems,  there  are  two 
practical  aspects  of  the  calculations  which  must  be  considered.  The  first 


concerns  the  choice  of  the  integration  framework  for  describing  the  molecular 
wavefunct ions  and  the  second  deals  with  the  choice  of  the  exchange  parameter, 
Qt,  in  different  regions  of  space. 

In  computations  with  heteronuclear  molecules,  there  are  several  free 
parameters  that  must  be  chosen:  the  ratio  of  sphere  radii  for  the  atomic 
sphe  res  of  integration  at  a  given  internuclear  separation,  the  degree  of 
sphere  overlap,  and  the  value  of  the  exchange  parameter  in  the  atomic  spheres 
and  the  intersphere  region. 

It  has  been  found  that  changing  the  ratio  of  the  sphere  radii  for  the  two 
atoms  in  a  heteronuclear  diatomic  molecule  introduces  changes  in  the  total 
energy  that  can  be  large  on  a  chemical  scale  (~  1  eV) .  A  choice  for  sphere 
radii  based  on  covalent  bonding  radii  does  not  necessarily  provide  a  good 
estimate  for  these  calculations.  The  value  of  the  exchange  parameter,  a,  and 
the  sphere  radii  and/or  sphere  overlap  is  normally  fixed  in  Xa  calculations 
for  crystals  where  the  geometry  is  fixed.  However,  to  develop  a  potential 
curve,  the  molecule  description  needs  to  change  substantially  as  the  inter¬ 
nuclear  separation  varies  and  the  changing  sphere  radii  include  varying  frac¬ 
tions  of  the  total  molecular  charge  (Reference  31).  Studies  made  at  UTRC  have 
shown  that  at  any  given  separation  the  total  energy  calculated  from  the  Xa 
model  is  a  minimum  at  the  radii  ratio  where  the  spherically  averaged  poten¬ 
tials  from  the  two  atomic  centers  is  equal  at  the  sphere  radius. 

''l'--,,’  ■  ''2‘'s2>  '2'‘> 

This  relationship  between  the  potential  match  at  the  sphere  boundary  and  the 
minimum  in  the  total  energy  appears  to  hold  exactly  for  "neutral"  atoms  and 
holds  well  for  ionic  molecular  constituents.  In  the  case  of  two  ionic 
species,  the  long  range  tail  of  the  potential  must  go  like  +2/R  from  one  ion 
and  -2/R  (in  Rydbergs)  for  the  other  ion  and  so  at  large  internuclear  separa¬ 
tions,  the  tails  of  the  potential  cannot  match  well.  However,  at  reasonable 
separations,  the  1/R  character  of  the  potential  does  not  invalidate  the  poten¬ 
tial  match  criterion  for  radii  selection.  This  match  for  the  atomic 
potentials  is  applied  to  the  self-consistent  potentials. 
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In  molecules  with  significant  charge  sharing  in  the  bonds,  the  radii  of 
the  atomic  spheres  is  frequently  increased  in  Xa  calculations  so  that  an  over¬ 
lap  region  appears  in  the  vicinity  of  the  bond  (Reference  32).  Studies  made 
at  UTRC  show  that  the  contribution  to  the  total  molecular  energy  from  the 
exchange  integral  shows  a  minimum  at  the  optimum  sphere  radius  or  sphere  over¬ 
lap.  This  provides  a  sensitive  criterion  for  selecting  these  parameters. 

The  values  of  the  exchange  parameters  in  the  spherical  integration  region 
around  each  atomic  center  are  frequently  set  at  the  atomic  values  both  for 


neutral  and  for 

ionic  molecular  const ituents - 

However,  for 

light  atoms,  the 

value  of  a  which 

best  reproduces  Hartree-Fock 

results  varies 

substant ial ly 

with  ionicity. 

In  argon. 

the  following  table 

compares ,  for 

the  neutral  atom 

and  the  positive 

ion,  the 

HF  energy  and  the  X 

a  energy  calcul 

ated  for  several 

values  of  a. 

a 

Xa  Energy 

HF  Energy 

Ar  ® 

.72177 

526.8176 

526.8173 

Ar^l/2 

.72177 

526.5857 

* 

.72213 

526.6007 

- 

Ar*  1 

.72177 

526.2447 

.72213 

526.2596 

- 

.72249 

526.2745 

526.2743 

The  optimum  value  of  a  changes  even  more  rapidly  in  the  fluorine  atom  goi 
from  0.73732  for  F®  to  .72991  for  F~^.  Since  the  total  energy  depends  li 
ly  on  a,  this  parameter  must  be  chosen  carefully. 

The  intersphere  exchange  coefficient  is  chosen  to  be  a  weighted  average 
of  the  atomic  exchange  parameters  from  the  two  constituents.  At  small  inter- 
nuclear  separations,  the  optimum  radius  for  an  atomic  sphere  frequently  places 
significant  amounts  of  charge  outside  that  atomic  sphere  -  charge  that  is 
still  strongly  associated  with  its  original  center  rather  than  being  trans¬ 
ferred  to  the  other  center  or  associated  with  the  molecular  binding  region. 

To  best  account  for  these  cases  the  weighting  coefficients  are  chosen  to 
reflect  the  origin  of  the  charge  in  the  intersphere  (or  outersphere 
region) , 
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where  (^Sj^  -  Q£°)  is  the  charge  lost  from  sphere  i  relative  to  its  atomic 

value  (or  ionic  value)  and  is  the  atomic  exchange  parameter  for  sphere 

i.  This  value  for  a.  ,  is  calculated  dynamically  -  it  is  updated 

intersphere  ■’ 

after  each  iteration  in  the  self-consistent  calculation. 

While  for  heavy  atoms,  these  changes  in  the  exchange  parameter  would  be 
small,  the  a's  for  small  atoms  vary  rapidly  with  z  (and  with  ionicity).  The 
correct  choice  of  the  exchange  parameters  influences  not  only  the  total  energy 
calculated  for  the  molecule  but  also  in  some  cases  affects  the  distribution  of 
charge  between  the  atomic  spheres  and  the  intersphere  region. 


2.2  TRANSITION  PROBABILITIES 


The  electronic  and  vibrational-rotational  wavefunct ions  of  a  pair  of 
states  can  be  used  to  calculate  transition  probabilities.  If  two  molecular 
states  are  separated  in  energy  by  an  amount  ~  hcv  (h  =  Planck's  constant 

c  =  velocity  of  light,  v  =  frequency  in  wave  numbers),  the  semi-classical 
theory  of  radiation  (References  33  and  34)  yields  for  the  probability  of  a 
spontaneous  transition  from  an  upper  state  n  to  a  lower  state  m 


^  nm  ■ 


_4 

3 


Snm 

gn 


(26) 


Here  A  is  the  Einstein  coefficient  for  spontaneous  transition  from  level  n 
nm 

m,  gj,  is  the  total  degeneracy  factor  for  the  upper  state 


gn  =  (2  -  8£r.A'){2S'+  ))(2J'  +  I) 


(27) 


and  S  is  the  total  strength  of  a  component  line  in  a  specific  state  of 
nm 

polarization  and  propagated  in  a  fixed  direction.  A  related  quantity  is  the 
mean  radiative  lifetime  of  state  n  defined  by 
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the  summation  being  over  all  lower  levels  which  offer  allowed  connections. 
The  intensity  of  the  emitted  radiation  is 


Inm"  ^^nm^n^nm 


(29) 


where  N^  is  the  number  density  in  the  upper  state  n.  This  analysis 
assumes  that  all  degenerate  states  at  the  same  level  n  are  equally  populated, 
which  will  be  true  for  isotropic  excitation.  The  total  line  strength  S  can 
be  written  as  the  square  of  the  transition  moment  summed  over  all  degenerate 
components  of  the  molecular  states  n  and  m: 


where  j  and  i  refer  to  all  quantum  numbers  associated  collectively  with 
upper  and  lower  electronic  states,  respectively. 

In  the  Born-Oppenheimer  approximation,  assuming  the  separability  of  elec¬ 
tronic  and  nuclear  motion,  the  wavefunction  for  a  diatomic  molecule  can  be 


written  as 


'^vjmA”  '^el  '/'v  (R)  X  > ) 


where  ib*’,  (r,  R)  is  an  electronic  wavefunction  for  state  i  at  fixed  inter- 
nuclear  separation  R,  iJ\,(R)  is  a  vibrational  wvefunction  for  level  v  and 
ipjn  <l>)  refers  to  the  rotational  state  specified  by  electronic  angular 

momentum  A,  total  angular  momentum  J  and  magnetic  quantum  number  M.  The 
representation  is  in  a  coordinate  system  related  to  a  space-fixed  system  by 
the  Eulerian  angles  (9,  x,  <|)) .  The  transition  moment  can  be  written, 

using  the  wavefunction  given  by  Equation  (31),  as 
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The  subscripts  e,  v  and  r  refer  to  the  electronic,  vibrational  and  rota- 
tional  wavefunct ions  and  M  and  M  are  the  electronic  and  nuclear  electric 
dipole  moments,  respectively.  Integration  over  the  electronic  wavefunction, 
in  the  Born-Oppenheimer  approximation,  causes  the  contribution  of  the  nuclear 
moment  M''  to  vanish  for  i  ^  j.  The  electronic  dipole  moment  can  be  written 
(References  34  and  35)  in  the  form 
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where  the  primed  coordinates  refer  to  the  space  fixed  system,  the  coordin¬ 
ates  r|^  refer  to  a  molecule-fixed  system  and  ^(0,  x»  is  a  group  rotation 
tensor  whose  elements  are  the  direction  cosines  related  to  the  Eulerian  rota¬ 
tion  angles  (9,  Xj  '^)  ■  Using  bracket  notation,  Equations  (32)  and  (33)  can  be 
combined  to  yield  for  the  transition  moment 


jv'J'A'M' 


The  matrix  elements  <J'  A'  M' |  (9,  x»  A"  M">  determine  the  group  se¬ 

lection  rules  for  an  allowed  transition  and  have  been  evaluated  for  many  types 
of  transitions  (References  36-38).  Summing  Equation  (34)  over  the  degenerate 
magnetic  quantum  numbers  M'  and  M"  we  have  from  Equation  (30) 
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where  j" /["  is  the  Honl-London  factor  (References  39  and  40) 


is  the  band  strength  for  the  transition.  Combining  Equations  (27),  (29)  and 
(35),  we  have  for  the  intensity  of  a  single  emitting  line  from  upper  level 
n: 
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where  Nj"  is  the  number  density  in  the  upper  rotational  state  J'  and 
4,  =  +  1)  is  the  electronic  degeneracy.  Taking  an 

average  value  of  v"J"  whole  band.  Equation  (37)  can  be  summed  to 

yield  the  total  intensity  in  the  (v',  v")  band: 


'  nv  J  _  4  /  Pmv" 
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where  N^,  "  ji  N ^ ,  is  the  total  number  density  in  the  upper  vibrational  level 
v'  and  where  we  make  use  of  the  group  summation  property 

,‘'-^l  =  (2J'+|)  (39) 


Comparing  Equations  (29)  and  (38),  we  have  for  the  Einstein  spontaneous 
transition  coefficient  of  the  band  (v',  v") 


r..nv'-|5  nv' 
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Similarly,  the  lifetime  of  an  upper  vibrational  level  v'  of  state  n  can  be 
written 


"  m«n  ^ 

where  the  summation  runs  over  all  v"  for  each  lower  state  m.  Equation  (40) 
can  be  cast  in  the  computational  form 
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where  Ae”  and  p”  are  in  atomic  units.  It  is  also  often  convenient  to 
m  V  m  V 

relate  the  transition  probability  to  the  number  of  dispersion  electrons  needed 
to  explain  the  emission  strength  classically.  This  number,  the  f-number  or 
oscillator  strength  for  emission,  is  given  by 
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The  inverse  process  of  absorption  is  related  to  the  above  development 
through  the  Einstein  B  coefficient.  Corresponding  to  Equation  (29),  we  have 
for  a  single  line  in  absorption 
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where  K(  v)  is  the  absorption  coefficient  of  a  beam  of  photons  of  frequency 


V  and 
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is  the  Einstein  absorption  coefficient  for  a  single  line.  Summing  over  all 
lines  in  the  band  (v",  v'),  assuming  an  average  band  frequency,  we  obtain 
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where  N^„  =  j„  Nj„  is  the  total  number  density  in  the  lower  vibrational  state 
v".  Corresponding  to  Equations  (42)  and  (43)  we  can  define  an  f-number  or 
oscillator  strength  for  absorption  as 
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In  computational  form.  Equation  (47)  becomes 
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where  AE  ”  and  p”  are  in  atomic  units.  Combining  Equations  (40)  and 
m  V  in  V 

(43)  and  comparing  with  Equation  (47),  we  see  that  the  absorption  and  emission 
f-numbers  are  related  by 
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Some  caution  must  be  observed  in  the  use  of  f-nurabers  given  either  by  Equa¬ 
tion  (43)  or  (47)  since  both  band  f-numbers  and  system  f-numbers  are  defined 
in  the  literature.  The  confusion  arises  from  the  several  possible  band  aver¬ 
aging  schemes  that  can  be  identified. 

An  integrated  absorption  coefficient  (density  corrected)  can  be  defined 
from  Equation  (46)  as 
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where  Che  exponential  factor  corrects  for  stimulated  emission.  Equation  50 
can  be  written  in  terms  of  the  absorption  f-number  as 
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Using  h  c/k  =  1.43880  cm-K* ,  we  obtain  a  computational  formula  for  the 

-2.  -1 

integrated  absorption  coefficient  as  S n  i  (cm  "atm  )  = 
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The  total  integrated  absorption  is  found  from 

SjOTAL^Z  ^  Sv",v' 
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(53) 


where,  under  normal  temperature  conditions,  only  the  first  few  fundamentals 
and  overtones  contribute  to  the  summations. 

The  developments  given  above  are  rigorous  for  band  systems  where  an  aver 
age  band  frequency  can  be  meaningfully  defined.  Further  approximations,  how¬ 
ever,  are  often  made.  For  example,  the  electronic  component  of  the  dipole 
transition  moment  can  be  defined  as 


This  quantity 
can  sometimes 
factored  form 


is  often  a 
be  chosen, 
as 


slowly  varying  function  of  R  and  an  average  value 
Equation  (36)  can  then  be  written  approximately  in 
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where  q  ,  „  the  square  of  the  vibrational  overlap  integral,  is  called  the 


Franck-Condon  factor.  ^ evaluated  at  some  mean  value  of  the  internuclear 
separation  R.  In  addition,  it  is  sometimes  possible  to  account  for  a  weak  R- 


©  .  . 
dependence  in  M  by  a  Taylor  series  expansion  of  this  quantity  about  some 

reference  value,  fi  „  usually  referred  to  the  (0,  0)  band.  We  have 

a  Py 
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Substituting  into  Equation  (56)  and  integrating  yields 
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where 
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is  the  R-centroid  for  the  transition  and 
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is  the  R^ceatroid.  Note  that  this  last  term  differs  (to  second  order)  from 
the  square  of  the  R-centroid.  An  alternate  procedure  can  be  developed  by 
evaluating  Equation  (54)  at  each  R-centroid,  Ryi^n.  Then 
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Equation  (60)  assumes  that  the  vibrational  wavef unction  product  \p^,  4)^„, 
behaves  like  a  delta  function  upon  integration, 


(61) 


The  range  of  validity  of  Equation  (60)  is  therefore  questionable,  particu¬ 
larly  for  band  systems  with  bad  overlap  conditions  such  as  oxygen  Schumann- 
Runge.  The  range  of  validity  of  the  R-centroid  approximation  has  been  exam¬ 
ined  by  Frazer  (Reference  41). 


The  final  step  in  calculating  transition  probabilities  is  the  determina- 


of  internuclear  separations,  R,  reached  in  the  vibrational  levels  to  be  con¬ 
sidered.  This  can  be  expressed  in  terras  of  the  expansion  of  Equation  (4) 
as 
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where  c  and  c^  are  coefficients  for  and  respectively. 

An  analysis  similar  to  that  yielding  Equation  (10)  and  (11)  gives 
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The  spatial  integral  in  Equation  (63)  reduces  to  one-electron  integrals 
equivalent  to  overlap  integrals,  and  the  evaluation  of  Equation  (63)  can  be 
carried  out  by  the  same  computer  programs  used  for  Equation  (11).  Programs 
for  evaluating  Equation  (62)  have  been  developed  at  UTRC  and  exam¬ 

ples  of  their  application  have  appeared  in  the  literature  (Reference  8). 

For  perturbed  electronic  systems,  the  transition  dipole  moment  will  have 
a  strong  R-dependence  and  R-centroid  or  other  apporximat ions  will  be  invalid. 
A  direct  evaluation  of  Equation  (36)  would  therefore  be  required  using  the 
fully-coupled  system  of  electronic  and  vibrational  wavefunct ions  to  properly 
account  for  the  source  of  the  band  perturbations. 


SECTION  3 


DISCUSSION  OF  RELATIVISTIC  METHODS 

For  heavy  atoms  (Z  >  30),  and  molecular  systems  built  from  heavy  atoms, 

relativistic  effects  become  increasingly  important  and  should  be  taken  into 

account  in  the  calculation  of  the  radial  wavefunctions .  The  implementation  of 

relativistic  effects  into  atomic  and  molecular  computer  codes  is  only  fairly 

recent  owing  to  the  increased  complexities  introduced  in  the  self-consistent 

field  (SCF)  procedure  and  the  greatly  increased  computer  time  required  for 

such  calculations.  Compared  with  the  non-relat ivist ic  case,  the  Dirac- 

Hartree-Fock  (DHF)  method  requires  that  two  radial  functions,  G  corre- 

ni  j 

sponding  to  the  large  component  and  F  .  corresponding  to  the  small  component 

J  1 

must  be  calculated  for  each  of  the  two  possible  j  values.  Thus,  the  numerical 
work  of  a  DHF  relativistic  treatment  is  increased  by  nearly  a  factor  of  four 
over  the  nonrelativistic  case,  exclusive  of  increased  complexities  in  evalua¬ 
tion  of  the  terms  of  the  Hamiltonian.  In  view  of  this,  methods  that  have  been 
developed  to  date  for  molecular  systems  have  involved  the  use  of  model  poten¬ 
tials  to  represent  relativistic  effects. 

In  the  calculation  of  the  internal  energy  of  a  molecular  system  comprised 
of  n  electrons  and  N  nuclei,  and  considering  only  electrostatic  interactions 
between  the  particles,  we  have  for  the  total  Hamiltonian 
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where  m  ,  m  ,  M  ,  are  the  masses  of  the  electron,  atom  a  and  combined  system 
e  a  T 

mass,  respectively.  Now  since  the  ratios  m  /m  and  m  /M,_  are  both  small, 


(2x10  -  5x10  )  we  can  effect  a  separation  of  the  electronic  and  nuclear 

coordinates  treating  the  total  wavefunction  as  a  product  of  a  nuclear  and  an 
electronic  part.  We  have 
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where  electronic  wavefunction  parametric  in  the  nuclear 

coordinates  as  given  in  Equation  (66)  and  nuclear  motion  wavefunc- 

tions  which  satisfy  (neglecting  terms  of  the  order  of  m^/m^) 
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The  cross  term  in  can  be  eliminated  by  a  proper  change  of  varia¬ 

bles  and  Equation  (67)  then  reduces  to  a  3N-3  dimensional  Schrodinger 
equation. 

For  most  systems,  where  the  velocity  of  motion  of  the  nuclei  is  slow 
relative  to  the  electron  velocity,  this  decoupling  of  electronic  and  nuclear 
motion  is  valid  and  is  referred  to  as  the  adiabatic  approximation.  Equation 
(66)  thus  defines  an  electronic  eigenstate  >  parametric  in  the 

nuclear  coordinates,  and  a  corresponding  eigenvalue  Ej^(Rjj)  which  is  taken  to 
represent  the  potential  energy  curve  or  surface  corresponding  to  state  k. 

In  the  usual  ab  initio  method  for  calculating  the  electronic  properties 
of  a  molecular  system,  one  starts  from  a  zero-order  Hamiltonian  that  is  exact 
except  for  relativistic  and  magnetic  effects,  and  which  involves  the  evalua¬ 
tion  of  electronic  energies  and  other  relevant  quantities  for  wavefunct ions 
that  are  properly  antisymmetrized  in  the  coordinates  of  all  the  electrons. 

For  a  system  containing  n  electrons  and  M  nuclei,  the  zero-order  Hamiltonian 
depends  parametrically  on  the  nuclear  positions  and  is  of  the  form 


1=1  1  =  1 
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where  z.  and  R.  are  the  charge  and  position  of  nucleus  i,  r.  is  the  position 

2  .  .  ^  . 
of  electron  j,  and  is  the  Laplacian  operator  for  electron  j.  All  quanti¬ 

ties  are  in  atomic  units,  i.e.  lengths  in  bohrs,  energies  in  hartrees  (1 
hartree  =  2  Rydbergs). 

In  addition  to  the  electrostatic  contribut ion,®^^ ,  the  complete  Hamilton¬ 
ian  should  contain  additional  terms  which  correct  for  magnetic  interactions 
and  relat ivis it ic  effects.  These  correction  terms  may  be  of  importance  in 
several  applications.  These  include: 

(1)  calculation  of  the  probability  of  making  a  transition  from  one  quan¬ 
tum  state  to  another  in  high-momentum  collisions  such  as  those  that 
can  occur  in  hot  atom  or  heavy  atom  chemical  dynamics  experi¬ 
ments; 

(2)  determination  of  the  interaction  energy  in  heavy  nuclei  systems  such 
as  Cs2  and  UO'*' ,  which  exhibit  open-shell  structure  on  both  nuclei  at 
infinite  internuclear  separations; 

(3)  calculation  of  the  intermolecular  forces  between  free  radicals, 
electronically  excited  states  of  molecules  with  open-shell  struc¬ 
ture,  and  long  molecular  conformations  of  possible  biological 
interest . 

3.1  BREIT-PAULI  HAMILTONIAN 

The  relativistic  correction  terms  to  the  usual  electrostatic  Hamiltonian 
have  been  derived  through  order  a  ,  where  o  is  the  fine  structure  constant, 
and  are  often  referred  to  as  the  Breit-Pauli  (Reference  42)  Hamiltonian  terms. 
This  Hamiltonian  has  been  derived  by  Bethe  and  Salpeter  (Reference  43)  for  a 
two-electron  system  and  has  been  generalized  to  the  many-electron  system  by 
Hirschfelder,  et  al  (Reference  44)  and  Itoh  (Reference  45).  In  the  absence  of 
external  electric  or  magnetic  fields  we  can  represent  these  correction  terms 
as  follows.  Let  s.  and  p.  =  -r  A .  denote  the  operators  for  the  spin  and  linear 


moment  of  electron  j,  respectively.  Then  the  generalized  Breit-Pauli  Hamil 
tonian,  correct  to  terms  of  0( a^/M) ,  can  be  written  as: 


*^BP  *^6  ‘"lL  ‘^SS  ‘"ls  *^0 
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where  is  given  by  Equation  (68)  and  the  correction  terms  can  be 
expressed  as  follows: 
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The  first  correctioa  term, ^  ,  represents  the  magnetic  orbit-orbit  coupling 

Li  Li 


terms  of  the  electrons  arising  from  the  interaction  of  the  magnetic  fields 

created  by  their  motion.  The  second  term,^  ,  gives  the  spin-spin  magnetic 

So 


coupling  terms  which  are  often  quite  appreciable.  For  r  =  0,  only  the 

J  ^ 


delta-function  contribution  survives  which  represents  the  Fermi-contact  spin 


interaction.  The  third  term,»^  is  usually  the  largest  in' magnitude  and 

LiS 


represents  the  spin-orbit  interaction  between  the  spin  and  magnetic  moment  of 
each  electron  and  the  spin-other  orbit  interaction,  which  represents  the 
coupling  of  the  spin  of  one  electron  with  the  magnetic  moment  of  a  different 
electron.  The  term.^p,  corrects  for  variation  of  the  electron  mass  with 
velocity  and  the  term,  represents  electron  spin  terms  identified  by 

Dirac  which  appear  to  have  no  classical  analogue. 


Aside  from  the  spin-orbit  term,^j^g,  usually  only  the  last  term,»^p. 


(often  called  the  Darwin  correction  term)  and»^p,  the  mass-velocity  term, 
are  retained  in  the  Hamiltonian,  yielding  the  so-called  Pauli  approximation 
(Reference  43). 

The  eigenfunctions  of  the  Hamiltonian  represented  by  Equation  (69)  are 
four-component  Dirac  spinors  which  may  be  expressed  as: 


(75) 


r  f 
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where  )^(  9,  <ti)  are  products  of  spherical  harmonics  and  Pauli  spinors  and 
P^j^(r),  Qj^j^(r)  represent,  respectively,  the  large  and  small  components  of  the 
radial  wave  equation.  The  exact  solution  of  the  Breit-Pauli  Hamiltonian  has 
only  been  given  for  one-  and  two-electron  atomic  systems  (Reference  43)  owing 
to  the  complexity  of  the  operators  for  the  general  n-electron  case.  For  a 
molecular  system,  Kolos  and  Wolniewicz  (Reference  46)  have  calculated  the 
relativistic  corrections  to  H2  using  Equation  (69)  to  (^(  a^) .  No  heavier  mole¬ 
cular  systems  have  been  treated  using  the  full  Breit-Pauli  Hamiltonian. 
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3.2  APPROXIMATE  TREATMENTS 


Although  the  Breit-Pauli  Hamiltonian  given  in  Equation  (69)  can  formally 
be  employed  in  a  molecular  system,  both  the  multiplicity  of  terms  and  the 
difficulty  of  evaluation  of  the  resultant  molecular  integrals  has  precluded 
its  general  use  to  date.  For  atomic  systems,  various  approximate  methods  of 
solution,  within  a  Hartree-Fock  or  multiconfiguration  Hartree-Fock  framework, 
have  been  proposed  for  atoms  (References  47-51).  In  most  of  these  methods,  a 
restricted  Hamiltonian  which  includes  only  the  one-electron  Dirac  terms  is 
usually  employed.  The  contributions  of  the  Breit  operators  for  spin-magnetic 
interactions  and  velocity  retardation  are  then  calculated  as  first-order  per¬ 
turbations  using  the  zeroth-order  Dirac  relativistic  wavefunct ions . 

An  even  more  approximate  method  for  incorporating  the  major  relativistic 
effects  has  been  proposed  by  Cowan  and  Griffin  (Reference  52).  In  this 
method,  the  mass-velocity  and  Darwin  Ud)  terms,  written  in  terms  of 

the  Pauli  equation  for  one-electron  atoms,  are  simply  added  to  the  usual  non- 
relativistic  Hamiltonian  operator.  In  addition,  the  spin-orbit  terms, 
are  omitted,  thereby  reducing  the  system  of  equations  to  a  single  form  repre¬ 
senting  the  description  of  the  major  component  wavefunction,  evaluated 

at  the  center-of-gravity  of  the  spin-orbit  states.  The  rationale  for  this 
approximation  lies  in  the  observation  that  detailed  atomic  calculations  using 
the  complete  DHF  method  have  indicated  that,  even  for  an  atom  as  heavy  as 
uranium,  less  than  1  percent  of  the  total  charge  is  described  by  the  small 
component  radial  wavefunct ions . 

The  resulting  equations  have  the  form; 


where  the  mass-velocity  and  Darwin  terms  take  the  form 


(77) 
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and  a  -  1/137.036  is  the  fine  structure  constant.  The  spin-orbit  term  is 
omitted  in  Equation  (76)  and  thus  these  equations  represent  center  of  gravity 
radial  functions  averaged  over  the  two  possible  total  angular  momentum  quantum 
numbers.  Equation  (76)  represents  (apart  from  the  neglect  of  spin-orbit 
effects)  the  relativistic  corrections  to  first  order  in  a^.  A  more  accurate 
analysis  of  heavy  atom  energy  levels  and  spectra  is  available  through  the  use 
of  the  radial  functions,  found  from  Equation  (76),  and  a  first-order 

perturbation  calculation.  Cowan  and  Griffin  (Reference  52)  have  illustrated 
the  utility  and  accuracy  of  such  an  approach. 

Recently,  Wood  and  Boring  (Reference  53)  have  adapted  this  approximate 
relativistic  method  to  the  local  exchange  problem  and  have  implemented  the 
solution  of  Equation  (76)  within  the  context  of  the  multiple  scattering  Xo 
method  (Reference  22).  The  central  field  Hamiltonian  is  modified  to  include 
mass-velocity  and  Darwin  terras,  given  by  Equations  (77)  and  (78),  in  the 
sphere  surrounding  each  atomic  center.  The  intersphere  region  in  the  multiple 
scattering  approach  (constant  potential  region)  is  treated  nonrelativistically 
since  charge  in  this  region  is  far  from  a  nucleus  and  is  screened  by  the 
charge  concentrated  around  the  atomic  centers.  The  matching  conditions  for 
continuity  of  the  wavefunction  at  the  sphere  boundaries  permits  any  necessary 
charge  transfer  between  the  relativistic  intra-atomic  regions  and  the  nonrela- 
tivistic  interatomic  constant  potential  regions.  For  an  atom,  the  Wood-Boring 
treatment  reduces  to  the  Dirac-Slater  local  exchange  method,  but  with  the 
neglect  of  spin-orbit  terras. 
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The  implementation  of  Equation  (76)  into  existing  nonrelativistic  multi¬ 
ple  scattering  molecular  codes  is  facilitated  by  a  change  in  the  dependent 
variable,  eliminate  the  first  derivative  of  the  wavefunction, 

illustrated  in  Equation  (78).  The  usual  Numerov  method  of  solution  can  then 
be  applied  to  the  central  field  problem;  the  only  new  requirement  being  the 
numerical  tabulation  of  the  first  and  second  derivatives  of  the  potential  at 
each  grid  point  in  the  integrations.  These  derivatives  are  computed  only  once 
for  each  complete  SCF  cycle  and  thus  the  total  required  computer  time  for  a 
typical  problem  is  not  significantly  increased  as  compared  with  a  nonrelativ¬ 
istic  calculation.  A  complete  self-consistent  program  incorporating  this 
method  has  been  developed  at  UTRC.  Our  code  has  been  tested  by  repeating 
calculations  for  the  U  and  Pu  atoms  (Reference  52),  where  we  find  excellent 
agreement  with  the  more  exact,  but  cumbersome,  Dirac-Slater  calculations. 
Results  for  molecular  calculations  have  recently  been  reported  by  Boring  and 
Wood  for  UF(,  and  U02^*  (References  54  and  55).  These  calculations  were  car¬ 
ried  out  to  illustrate  the  shifts  in  the  valence  levels  for  such  systems  re¬ 
sulting  from  relativistic  effects.  The  total  energy  was  not  of  principal 
concern. 

We  have  recently  reported  (Reference  56)  the  first  all-electron  calcula¬ 
tion  of  the  potential  energy  curves  for  a  molecule  (Hg^^)  built  from  atoms 
which  exhibit  significant  relativistic  effects.  This  study  illustrated  that 
reliable  total  energies  are  obtainable  through  a  relativistic  multiple  scat¬ 
tering  density  functional  treatment,  provided  care  is  taken  to  optimize  poten¬ 
tial  match  and  overlap  criteria  for  such  systems.  This  study  formed  the  basis 
of  the  computational  scheme  that  we  have  employed  here  for  the  uranium/oxygen 
system. 

3.3  EFFECTIVE  CORE  MODELS 

It  is  well  known  from  chemical  experience  that  the  outermost  valence 
electrons  contribute  most  to  determining  the  chemical  properties,  especially 
spectroscopic  properties,  of  molecules.  The  core  electrons  remain  essentially 
unchanged  from  their  atomic  form  except  for  internuclear  separations  of  the 
order  of  the  charge  radii  of  the  outer  core  region  or  less,  wherein  core  po¬ 
larization  effects  may  become  important.  Since  the  computational  time 


required  for  ab  initio  calculations  of  electronic  structure  goes  up  at  least 
quadratically  with  the  number  of  electrons  in  the  system,  there  have  been  many 
attempts  to  replace  the  more  tightly  bound  core  electrons  with  simple  one- 
electron  effective  potentials  (References  57-65).  Concurrent  with  elimination 
of  an  explicit  treatment  of  the  core  electrons,  a  transformation  of  the 
valence  orbital  basis  is  required  to  insure  that  the  lowest  valence  orbital  of 
each  sjnranetry  has  a  nodeless  radial  form,  since  it  is  well  known  that  the 
lowest  energy  eigenfunction  for  a  local  potential  must  be  nodeless  (Reference 
66). 

Typical  of  the  several  effective  core  models  that  have  been  reported  is 
that  due  to  Kahn,  et  al  (Reference  64)  whereby  an  effective  core  potential  is 
described  in  terras  of  angularly  dependent  projection  operators  as 

u“'' =  f  I  I  |ltn>[u,“''(r)-u,'°''(,)]<im| 

where  L  is  taken  at  least  as  large  as  the  highest  angular  momentum  orbital 
occupied  in  the  core.  The  terra  (r)  represents  the  effective  Coulomb  and 

exchange  potential  felt  by  the  valence  electrons.  The  second  term  essentially 
accounts  for  the  repulsive  potential  between  valence  and  core  electrons  for 
each  symmetry  I,  The  only  non-local  character  exhibited  by  a  potential  of  the 
form  of  Equation  (79)  arises  from  the  Z-dependence  which  can  be  cast  in  terms 
of  one-electron  integrals  between  the  core  and  valence  orbitals.  Explicit 
two-electron  terms  connecting  core  and  valence  orbitals  are  thus  avoided  which 
greatly  simplifies  the  calculation  of  matrix  elements  of  the  effective  Hamil¬ 
tonian.  The  potential  given  by  Equation  (79)  can  be  compared  with  the  gener¬ 
alized  Phillips-Kleinman  pseudo-potential  (Reference  59). 

(80) 

^core  .  ^  ^  yCO 

where  J^,  and  K(.  represent  the  core  orbital  Coulomb  and  exchange  operators 
CO 

and  V  is  a  complicated  non-local  operator  which  guarantees  core-valence 
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orthogonality.  Since  the  P-K  core  orbitals  must  simultaneously  be  eigenfunc- 

CO 

tions  of  both  the  core  and  valence  Hartree-Fock  Hamiltonians,  V  ,  in  general, 
contains  complicated  two-electron  terms  and  limits  the  usefulness  of  Equation 
(80)  over  a  full  ab  initio  treatment. 

The  prescription  of  Kahn  can  be  implemented  by  analytically  fitting  a 
nodeless  pseudo-orbital,  to  a  linear  combination  of  numerical  or  analytic 

Hartree-Fock  orbitals  determined  from  a  full  self-consistent  treatment  of  the 
core  electrons,  using,  for  example,  the  multiconfiguration  Hartree-Fock  code 
of  Froese-Fischer  (Reference  67).  The  components,  U  (r),  of  Equation 

«•  X* 

(79)  are  then  defined  implicitly  from  the  Schrodinger  equation 
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Equation  (81)  can  be  extended  to  relativistic  systems  in  several  ways. 
Kahn,  et  al  (Reference  68)  suggest  an  approximate  treatment  of  adding  only 
the  mass-velocity  and  Darwin  terras  to  the  usual  electrostatic  Hamiltonian  and 
to  determine  approximate  HF  orbitals  in  the  manner  prescribed  by  Cowan  and 
Griffin  (Reference  52).  Equation  (81)  is  then  used  to  determine  an  effective 
(r)  such  that  are  the  eigenvalues  of  the  CG  approximate  relativis¬ 
tic  solution  and  are  curve-fitted  to  the  CG-HF  orbitals.  In  this  treat¬ 

ment,  the  represent  approximate  solutions  to  the  major  component  wavefunc- 
tion,  P„jj,  determined  at  the  center-of-gravity  of  the  spin-orbit  states. 

Lee,  et  all  (Reference  65)  adopt  a  somewhat  more  complicated  treatment  in 
which  the  spin-orbit  operator  is  added  to  the  usual  electrostatic  Hamiltonian, 


in  addition  to  the  mass-velocity  and  Darwin  terms  retained  in  the  Cowan- 
Griffin  treatment.  The  large  component  eigenfunctions  of  a  full  Dirac- 
Hartree-Fock  treatment  of  the  atom,  as  given,  for  example,  by  Desclaux  (Refer¬ 
ence  69)  are  then  curve-fitted  in  a  manner  similar  to  the  Kahn  treatment  but 
include  the  additional  index  for  the  particular  spin-orbit  state,  Use 

of  these  eigenfunctions  in  a  molecular  system  fits  more  naturally  into  a  (J-J) 
coupling  scheme  whereas  the  determined  using  the  Kahn  method  are  more 
easily  represented  using  A-s  coupling. 

Although  these  effective  core  models  can  often  accurately  describe  an 
atomic  eigenvalue  sequence,  including  even  high-lying  electronically  excited 
states  (Reference  70),  there  are  inherent  difficulties  in  their  application  to 
molecular  environments,  where  the  maximum  angular  momentum  component  of  the 
valence  shell  orbitals  may  often  exceed  the  highest  £.-value  component  retained 
in  Equation  (82).  This  is  particularly  true  for  valence  orbitals  which  exhi¬ 
bit  strong  changes  from  atomic  form  through  hybridization  with  higher  angular 
momentum  orbitals  or  through  the  addition  of  more  compact  polarizaiton  terms. 
In  either  case,  since  the  relativistic  terms  are  now  all  buried  in  a  fixed 
rather  than  a  dynamic  relativistic  operator,  only  static  core  effects  are 
imposed  in  determining  the  shape  of  the  valence  molecular  orbitals.  Relativ¬ 
istic  effects  between  valence  electrons  and  shielding  effects  of  the  core  by 
the  valence  electrons  are  therefore  neglected  in  these  effective  core  treat¬ 
ments.  In  addition,  the  models  obviously  break  down  completely  when  the 
nuclei  are  brought  together  to  dimensions  such  that  core  overlap  and  polariza¬ 
tion  effects  become  significant.  Unfortunately,  calculations  to  date  seem  to 
indicate  that  such  effects  begin  to  set  in  for  internuclear  separations  of  the 
order  of  equilibrium  bond  lengths. 
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SECTION  4 


DISCUSSION  OF  RESULTS 

An  analysis  of  the  electronic  structure  of  UO  and  UO"*"  using  a  relativis¬ 
tic  formulation  has  previously  been  undertaken  by  UTRC  (Reference  76).  The 
possible  low-lying  molecular  states  are  shown  in  Tables  1  and  2.  Preliminary 
calculations  were  performed  for  several  states  of  UO  and  UO'*'  and  for  the 
ground  state  of  UO2,  UO2*  and  UO2**.  A  brief  summary  of  the  results  of  these 
molecular  calculations  follows. 

UO 

Electronic  structure  calculations  were  carried  out  for  this  system  using 
a  relativistic  density  functional  formalism.  Only  a  selected  group  of 
symmetries  was  studied.  Our  calculations  indicate  that  the  lowest  symmetry  of 
UO  is  derived  from  the  (A,S)  coupled  ^I  state  and  has  the  following  principal 
molecular  orbital  occupancy: 


*1  [lo^  2a  3a  4a^  In  2n  3n  U  H  5o] 
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The  1  3w  1  6  1  (ji  1  group  derives  from  the  5f  atomic  configuration  in  the  U 
atom  and  is  quartet  coupled.  We  have  found  that  the  ^I  state  of  UO  is  the 
ground  state  but  that  several  other  symmetries,  including  %  and  ^I,  are  low- 
lying.  Our  results  indicate  that  only  states  of  triplet  or  quintet 
multiplicity  are  bound  for  this  system.  An  examination  of  the  structure  of 
these  low-lying  states  of  UO  indicates  a  near  total  charge  transfer  to  U'*'^"^ 
for  short  (equilibrium)  internuclear  separations.  Thus  only  the  molecular 
states  in  the  lower  group  shown  in  Table  1  are  likely  to  be  bound.  This  would 
yield  42  bound  molecular  states  arising  from  u[^]  +  o[3p]  and  39  repulsive 
states.  Since  (J-J)  coupling  is  surely  a  better  approximation  for  UO,  these 
two  manifolds  of  states  will  be  optically  connected  and  many  pre-dissociation 
paths  of  the  type: 


UO**  +  hv  +  UO**  -*•  U  +  0 


(84) 
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are  possible.  Here  UO  is  a  vibrat ional ly  excited  low-mult iplet  state  of  UO 
** 

and  UO  is  a  dissociating  state.  The  predicted  optical  absorption  should  be 
strong  since  the  transfer  is  from  highly  ionic  states  to  neutral  valence 
states  of  UO.  Since  UO  has  a  large  dissociation  energy  (7.87  eV),  both  one 
photon  and  two  photon  solar  excitation  processes  are  possible. 

A  vibrational  analysis  of  the  0=5  ground  state  of  UO  was  carried  out 
using  a  Hulbert-Hirshfelder  (Reference  77)  fit  to  our  calculated  potential 
curves.  This  fit  yields  an  equilibrium  internuclear  distance  of  1.89  A  and  a 
fundamental  vibrational  constant  (o^)  of  859  cm“ The  spin-orbit  interac¬ 
tion  was  calculated  using  U'*’^  atomic  splitting  parameters.  No  explicit  two- 
center  effects  are  included.  Our  calculated  spectroscopic  data  are  compared 
in  Table  3  with  the  theoretical  work  of  Krauss  and  Stevens  (Reference  78),  the 
recent  rotational ly  resolved  experimental  studies  of  Heaven  and  Nicolai 
(Reference  79),  and  estimates  based  on  experimental  data  for  similar  systems. 
The  agreement  is  well  within  the  uncertainty  of  the  calculations  or  experi¬ 
mental  estimates.  The  theoretical  studies  predict  a  5  ground  state  for  UO 
whereas  the  experimental  data  of  Heaven  and  Nicolai  suggest  a  ^  ground 
state.  The  character  of  all  of  the  low-lying  multiplets  of  UO  is  similar 
however,  and  this  apparent  discrepency  should  not  affect  our  conclusions  about 
optical  or  LWIR  absorptions. 

An  analysis  of  the  LWIR  emission  from  UO  was  carried  out  based  on  the 
ground  electronic  state.  These  calculations  should  also  be  representative 
of  the  LWIR  emission  from  other  low-lying  electronic  states  since  they  exhibit 
similar  ionicities.  Our  calculated  f-numbers ,  including  fundamentals  and 
overtones,  were  reported  in  DNA-TR-82-159 .  These  data  (v'  >  v")  were  given 
for  the  lowest  30  vibrational  levels.  For  UO,  our  calculated  f-number  for  the 
1-0  transition  is  4.86  x  10~^  at  X  =  12.04  u. 

UO'*' 

Detailed  searches  of  several  symmetries  of  VO*  were  carried  out  to  deter¬ 
mine  the  ground  molecular  state  of  this  system.  Our  calculations  indicate 
that  the  lowest  sjrmmetry  of  UO"*”  is  derived  from  the  (A,  S)  coupled  **I  state 
and  has  the  follow'ng  principal  molecular  orbital  occupancy: 


[lo^  2o^  3cj^  4a^  lir^  Im*  }  3it  1(|(  16  I] 


(85) 


I  I 

Again  we  found  that  the  |  3n  1(^  16  group  is  quartet  coupled  in  the  ground 
state  but  a  second  manifold  of  states  for  the  UO'*'  system,  which  exhibit  doub¬ 
let  coupling  of  these  electrons,  lies  about  2  eV  above  the  ground  state.  An 
apparent  gap  in  the  density  of  states  for  UO'*’  is  found  between  these  two 
groups . 

Calculations  for  UO'*’  proved  to  be  much  more  complex  than  UO  owing  to  the 
presence  of  at  least  two  low-lying  dissociation  limits  of  U'*’  +  0.  As  indi¬ 
cated  in  Table  2,  the  ionic  U^'*’0“  structures  will  mix  in  all  multiplicities 
with  molecular  ion  states  arising  from  U^[^l]  +  0[^P],  U*  [^L  ]  +  o[^p],  and 
U'*'  [^K  ]  +  0  pP  ]. 

A  vibrational  analysis  of  the  Q  =  9/2  ground  state  of  UO'*'  was  carried  out 
using  a  Hulbert-Hirschfelder  (Reference  77)  fit  to  our  calculated  potential 
curves.  The  spin-orbit  splittings  were  derived  from  atomic  parameters  for  the 
O'*"  ^  ion.  This  fit  yields  an  equilibrium  internuclear  distance  of  1.84  A  and  a 
fundamental  vibrational  constant  of  890  cm" These  data  are  compared  in 
Table  3  with  other  calculated  estimates,  since  there  are  no  experimental  data 
available.  The  agreement  between  our  work  and  that  of  Krauss  and  Stevens 
(Reference  78)  is  less  satisfactory  than  in  the  case  of  UO  but  still  well 
within  the  uncertainty  of  the  several  calculations. 

A  perturbative  treatment  for  calculating  the  density  of  states  in  uranium 
molecules  is  available  through  the  use  of  ligand  field  theory.  The  basic 
concept  relies  on  the  assumption  that  the  structure  and  density  of  the  5f 
electrons  in  uranium  and  its  ions  remain  unchanged  in  a  molecular  environment. 
Our  calculated  spin-orbit  splitting  for  the  5f  electrons  in  U'*’ ^  +  U'*’ ^  is  given 
in  Table  4  where  a  comparison  with  earlier  studies  by  Wood  and  Boring 
(Reference  53)  is  given.  The  overall  agreement  is  excellent.  In  Table  5,  we 
present  the  calculated  excitation  energies  of  several  configurations  of  U'*’^. 
These  data  clearly  indicate  that  the  **1  component  of  5f  ^  is  low-lying  and  that 
the  lowest  doublet  manifold  lies  ~  2.2  eV  above  the  ground  state.  If  this 


/.  £. 


splitting  carries  over  to  the  UO  ion  without  much  change,  the  following  solar 
pumped  process  is  possible: 


UO'^[^l]  +  hv  -»•  UO'^*[^h] 


(86) 


e  +  UO'^*[^h]  -»■  U  [^l]  +  0[^P] 

L  j  l  J  L  J 

Dissociative-recombination  is  not  energetically  possible  from  the  ground 

+  «  •  «  'y 

state  of  UO  but  it  is  possible  from  the  excited  H  and  higher  states.  The 

optical  connections  between  the  manifolds  of  UO'*’  states  will  be  studied  in 

future  work. 

An  analysis  of  the  emission  characteristics  for  the  ground  state  of  UO"*” 
indicates  an  oscillator  strength  for  emission  (fj^o)  of  5.17  x  10~^  at  X  = 

11.3  U.  A  complete  analysis  of  our  calculated  LWIR  emission  for  UO'*'  was  given 
in  DNA  TR-82-159.  Our  calculated  LWIR  emission  for  UO"*"  is  typical  of  that  for 
a  highly  ionic  metal  oxide.  We  predict  strong  emission  from  the  fundamentals 
of  UO'*'  in  the  wavelength  region  11.3  ^  14  y.  Since  this  system  exhibits  weak 
anharmonicity ,  we  find  the  overtones  down  in  intensity  by  several  orders  of 
magnitude.  However,  the  first  excited  state  of  UO'*’  (**8)  lies  at  ''  2800  cm“^ 
in  our  calculations  with  a  predicted  electronic  oscillator  strength  of  ''  1  x 
10“^  for  the  **1  -  ‘*H  transition.  The  electronic  and  vibrational  manifolds  for 
UO^  are  highly  overlapped  above  the  second  vibrational  level  of  the  ground 
^lg/2  state.  Since  the  density  of  electronic  states  of  UO^  is  very  large 
above  ^2.2  eV ,  we  predict  that  strong  solar  pumping,  followed  by  intense 
radiation  in  the  region  0.6  <  X  <  11.3  y  should  occur  for  this  system.  This 
conclusion  is  similar  to  that  reached  by  Krauss  and  Stevens  (Reference  78) 
based  on  their  MCSCF  analysis  of  the  UO'*’  system.  Since  the  excited  electronic 
states  of  UO'*’  lying  in  the  region  of  strong  solar  flux  (  .4  -  .7  y)  exhibit 
shifted  equilibrium  internuclear  separation  from  that  of  the  ground  **I  state, 
we  predict  efficient  conversion  of  solar  photons  to  IR  photons  for  this 
system. 
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k  -  k  •  «  -i.'-  .V  .V  v*.  A 


U02^ 

An  analysis  of  U02^  has  been  carried  out  in  symmetry.  McGlynn  and 

Smith  (Reference  80)  have  proposed  a  ^4^  ground  state  for  this  ion  based  on 
simple  molecular  orbital  arguments.  Their  set  of  one-electron  orbital  ener¬ 
gies  is  based  on  a  maximum  overlap  criterion  that  is  empirical  in  character. 
More  modern  calculations  of  the  actinide  series  atoms  suggest  that  a  U  ^[o 
or  ^ ]2  structure  should  be  the  most  stable  configuration.  In  terms  of 

MO's,  the  lowest  predicted  electronic  state  (Dooj,)  would  be: 

log^  lou^  2og^  iTTu^  2au^  3ag^  2itu^  litg^  3au^  |  Hu  :  (88) 


2 


An  extensive  series  of  calculations  of  the  possible  low-lying  symmetries  of 
U02^  now  definitely  establish  this  to  be  the  ground  state,  but  also  predict  a 
very  low-lying  state. 

A  series  of  calculations  of  the  low-lying  electronic  states  of  U'*’** 

and  U'*'^  was  first  undertaken  to  determine  the  approximate  location  of  the 
electronic  states  of  the  central  ion  which  corresponded  to  f  d  transitions. 
These  results  are  given  in  Table  6  which  indicate  that  f  d  promotion  in 
these  ions  lies  at  327  ,  167  and  108  nm,  respectively,  for  U'*’^,  U'*'**  and  U'*’^. 
This  strong  absorption  (dipole-allowed),  corresponding  to  central  ion  promo¬ 
tion,  is  predicted  to  occur  only  for  wavelengths  shorter  than  ■“  300  nm  for 
U02^-  A  spin  orbit  analysis  of  these  calculated  levels  yields  the  splittings 
shown  in  Figure  1.  The  ‘®5/2u  predicted  ground  state  multi- 

plet . 

In  order  to  establish  that  f  -*■  d  transitions  in  U02^  lie  at  these  short 

wavelengths,  a  series  of  molecular  calculations  of  the  low-lying  electron 

states  was  undertaken.  The  promotion  of  an  electron  out  of  the  H^  MO  of 

002^,  which  is  nearly  purely  0  in  composition,  gives  rise  to  the  one-electron 

spectra  shown  in  Table  7.  The  u  +  u  intra-f-shel 1  transitions  are  dipole- 

forbidden  and  correspond  to  the  several  weak  absorption  bands  observed  in 

solutions  of  U02^  (Reference  81).  The  lowest  dipole-allowed  f  d  transitions 

on  the  central  uranium  atom  corresponds  to  the  (^f  )  +  ^A  (6d  )  excita- 

u  u  g  g 

tion  at  X  ^  300  nm.  This  transition  lies  beyond  the  region  of  efficient  solar 
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pumping  of  UO2''’ .  Hay  (Reference  82)  reports  a  similar  transition  in  his 
effective  potential  studies  of  U02^,  but  at  a  somewhat  shorter  wavelength 
C''  250  nm)  . 

A  separate  molecular  spectra  for  UO2''’  arises  from  charge  transfer  states 
formed  by  promotion  of  an  electron  from  (mainly)  ligand  MO's  to  a  central 
uranium  atom  MO.  Preliminary  calculations  (Reference  83)  indicated  that 
several  of  these  states  lie  in  the  300-400  nm  region,  and  would  therefore  be 
efficiently  solar  pumped.  A  careful  re-analysis  of  the  ground  ^4^  state 
indicated  a  false  convergence  onto  local  minima  in  the  orbital  parameter 
space.  This  produced  a  ground  state  potential  curve  with  an  internuclear 
separation  that  was  too  large,  yielding  excitation  energies  that  were 
correspondingly  too  small.  Careful  convergence  studies  and  the  employment  of 
a  new  integration  scheme  now  gives  an  equilibrium  U-0  bond  length  in  U02^  of 
1.73  A.  This  value,  coupled  with  a  fundamental  Vj  vibrational  frequency  of 
880  cm~^,  is  now  in  good  agreement  with  Badger's  rule  (Reference  84) 
correlating  bond  lengths  and  force  constants  for  actinide  salts  and  0x0- 
ions : 

R(U-O)  *  1.08  k  +  1.17  (89) 

Using  these  more  reliable  integration  methods,  the  low-lying  charge 

transfer  states  of  U02^  were  re-examined;  their  term  levels  are  given  in 

Table  8.  These  charge  transfer  states,  some  of  which  would  exhibit  very 

strong  absorption,  begin  at  ^  47  ,000  cra~^  [**n  ],  with  the  first  strong  dipole 

.  4  3  2 

allowed  transitions  corresponding  to  the  excitation  lirg  li|>u  lug  1  <t>u  . 

These  transitions  lie  at  53,000  57,000  cm  which  correspond  to  absorption 

wavelengths  in  the  VUV ,  well  outside  of  the  region  for  efficient  solar  pump¬ 
ing.  A  composite  potential  energy  curve  for  U02^,  including  both  central  ura¬ 
nium  atom  excitation  states  and  the  low-lying  ligand  to  central  atom  charge 
transfer  states,  is  given  in  Figure  2.  Our  calculated  assignments  of  the 
f  f ,  f  ->■  d  and  charge  transfer  excitations  among  the  low-lying  states  are 
given  in  Table  9. 

This  oxo-ion  exhibits  two  different  characteristic  equilibrium  U-0 
separations.  They  correspond  to  U*^  (0”^)  configurations  with  a  short  bond 


and  to  U'*’ **  (0~  2  charge  transfer  states  which  exhibit  a  weakened  U-0  bond 

and  correspondingly  longer  equilibrium  bond  length.  The  spin  orbit  analysis 
for  UO  2^  is  presently  being  carried  out  but  this  analysis  will  not  affect  our 
conclusion  that  solar  pumping  of  electronically  excited  states  of  UO  2'''  has  a 
low  efficiency. 

An  analysis  of  the  LWIR  emission  from  UO2''’  has  been  carried  out  based  on 
the  ground  electronic  state.  These  calculations  should  also  be  repre¬ 

sentative  of  the  LWIR  emission  from  other  low-lying  electronic  states  since 
they  all  exhibit  similar  ionicities.  The  spectroscopic  data  for  the  ^'^5^2u 
state,  derived  from  our  calcuated  ground  state  potential  curve,  are  shown  in 
Table  10.  We  find  a  small  anharmonic ity  in  the  IR  active  asymmetric  vibra¬ 
tion,  similar  to  that  found  in  the  highly  ionic  UO'*'  and  UO  molecules.  Our 
calculated  f-numbers,  including  fundamentals  and  overtones  are  given  in  Table 
11  for  the  lowest  30  vibrational  levels.  We  have  also  included  the  emission 
wavelengths  for  each  transition.  For  UO2'*',  our  calculated  f-number  for  the  1- 
0  transition  is  1.2  x  10“**  at  X  =  11.43  u.  This  can  be  compared  with  a  cal¬ 
culated  f-number  of  5.2  x  10”^  at  A  =  11.3  p  for  UO'*’.  Relatively  strong  LWIR 
is  then  predicted  for  002'*’.  In  Table  12,  we  present  the  calculated  integrated 
band  absorption  coefficients  as  a  function  of  temperature.  All  overtone 
contributions  have  been  included  in  these  band  absorption  coefficients. 


An  extensive  set  of  calculations  for  the  ground  state  of  UO^  was 
carried  out  in  symmetry.  The  ground  electronic  state  has  the  dominant 
molecular  orbital  configuration: 

[l  og^  1  ou^  2ag^  Ittu*^  2au^  3og^  2110*^  1  itg^  3ou^  1  ]  (90) 

We  find  an  equilibrium  internuclear  separation  of  1.81  A,  a  value  somewhat 
larger  than  that  corresponding  to  the  various  uranyl  oxo-ions.  At  present  it 
is  uncertain  whether  this  represents  a  deficiency  in  the  calculations  or 
whether  significant  back-bonding  of  electron  charge  to  the  central  uranium 
atom  is  occurring  in  the  gas  phase  UO2  species.  This  would  result  in  a 
lengthing  of  the  U-0  bond  and  a  corresponding  decrease  in  the  U-0  bond 
strength. 


wi  VJ  j-j-wsirj  trr-rj  -r.  w,-ir.  ^j-fjT.-wv'iT.'V’/^JirjTri^d  irvr^wjirirnwK.w  vnv^-ywrWuirsrr^ 


The  first  excited  electronic  state  of  UO2  corresponds  to  a  l(tiu  •*•  Situ 
electron  promotion.  This  state  lies  at  ■“  1.4  eV  but  we  predict  very  weak 
solar  pumping  since  the  transition  is  not  dipole  allowed.  The  first  strongly 
allowed  transition  corresponds  to  a  l^iu  -*■  l6g  electron  promotion  at  ''  3.4  eV. 
This  transition  lies  at  360  nm  putting  it  out  of  the  region  of  efficient 
solar  pumping.  The  LWIR  analysis  for  neutral  UO2  should  be  carried  out, 
however,  since  the  highly  ionic  nature  of  UO2  may  give  rise  to  a  strong 
absorption/emission  character  for  the  asymmetric  stretch  vibrational  mode. 
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SECTION  5 


RECOMMENDATIONS 

The  low-lying  molecular  states  of  UO'*’  arising  from  **I  (5f^7s^), 

^L(5 f ^7 s6d) ,  and  ^K(5f^7s6d)  of  are  given  in  Table  2.  Additional  low-lying 
states  of  UO'*’  arise  from  the  ^M(5f^6d7p)  state  of  U'*’  and  from  the  doubly 
ionized  ^L(5f^6d)  and  ^I(5f^7s)  and  triply  ionized  **I(^f^)  states  of  the 
uranium  ion.  Further  calculations  of  the  term  manifold  for  this  molecular  ion 
and  doubly  excited  states  of  U'*’  and  UO'*’  are  needed  in  order  to  analyze  the 
role  of  dielectronic  recombination  in  this  system.  In  our  preliminary 
studies,  a  separate  (higher  lying)  chemistry  was  found  for  U'*’  when  the  5f 
electrons  were  doublet  spin-coupled.  Verification  of  this  structure,  using  a 
relativistic  framework  for  this  ion,  is  an  area  for  future  research  activi- 
t  ies . 

Follow-on  studies  will  also  include  an  analysis  of  the  several  possible 
kinetic  pathways  that  have  been  suggested  by  our  work  to  date.  These  include 
photoexcitation  and  photodissociation  of  UO  as: 

U'^O"  +  hv  -*•  UO* 

*  (91) 

UO  +  hv  +  U  +  0 

where  UO*  represents  the  excited  neutral  valence  states  of  UO .  Other 
potentially  important  kinetic  studies  include  dissociative-recombination 
routes  from 

e  +  [UO^]*  -»■  U*+  0  (92) 

where  [uO'*’ ]  is  an  electronically  excited  state  of  UO'*'  which  exhibits  doublet 
coupling  for  the  5f  electrons  and  charge  transfer  processes  such  as 


In  addition,  the  dissociative-recombination  routes 


e  +  > 

e  +  UO^ 


U"^  +  0 
U  +  0^^ 

U'*'  +  0  +  0 
UO  + 

UO  +  0 


(94) 


are  uncertain  until  the  energetics  of  the  doubly  ionized  UO’*"*’  and  1102^''’ 
species  are  better  defined. 

In  addition  to  these  proposed  calculations  of  the  structure  and  radiative 
properties  of  uranium/oxygen  species,  an  analysis  of  the  importance  of 
dielectronic  recombination  processes  should  be  carried  out.  Dielectronic 
recombination  (DR)  is  the  process  by  which  electron  capture  from  the  continuum 
to  a  bound  state  is  facilitated  by  excitation  of  a  previously  bound  electron. 
The  process  can  be  represented  as: 


e  +  >  [uf* 

[Uf*  >  [uf  +  hv 


(95) 


where  [u  ]**  is  a  doubly  excited  continuum  state  of  the  neutral  system  and 

[u  ]*  is  a  singly-excited  Rydberg  or  Rydberg-like  structure.  A  corresponding 

process  is  possible  for  e  +  UO"''.  The  low-lying  doubly  excited  state  structure 

of  uranium  begins  at  about  4600  cm”^,  indicating  that  this  process  needs  to  be 

considered  for  T  >  3000  "K.  The  region  where  e  +  UO^  is  important  is  unknown 
e 

at  the  present  time. 


(AXIAL)  LIGAND  FIELD 
SPLITTING 


CALCULATED  SPLITTING 


Figure  1.  UO^  energy  levels  —  Sf  splitting. 
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Table  1.  Molecular  states  of  UO. 


Separated  Atom 
U  +  0 


♦  ^Pg(2p^) 

£•  -  0.000  eV 
(81  states) 


Molecular  States 


^^■(1).  ^r(i),  ^r‘*’(2),  ’ 

’£■^(2),  ^n(3).  ^n{3),  ’n(3),  ^a{3 
^A(3).  ^♦(3),  ^*(3),  ^♦O),  ^r(3) 

^r(3),  3h(3),  5h(3).  ^H(3),  ^1(3) 

^1(3),  hO),  V(3),  \(3),  ^A(2) 

^A(2),  ^M(l),  ^M(l),  ’m(I) 


Ky(5f^7s^6d)  +  ^Pg(2p^)  ^^•(2),  V(l),  ^£"(2),  ^£'^{1), 

E.  -  0.096  eV  ^n(3),  ^n(3),  ^n(3),  ^A(3),  ^A(3)  , 

(72  states) 


Separated 
U"*"  +  0~ 

\i5thsh  +  ■ 

(42  states) 

U'*"*'  +  0° 
5Ly(5f^6d)  +  ^ 


Table  1.  Molecular  states  of  UO  (Concluded) 

Atom  Molecular  States 


P^(2p5)  h-i2),  h-(2),  % 

5n(3),  3a(3),  5a(3),  3*0) ,  5*0)^ 
3r(3),  3r(3),  3h(3),  5h(3),  ^1(2), 
5i(2),  3k(1),  5k(1) 

lg(2p6) 


(9  states) 


5e-(1),  5n(i),  5^(1),  5*(i)^  5r(i) 
5h(1),  5i(i),  5k(1),  5/y(i) 


Table  2.  Molecular  states  of  UO"*". 


Separated  Atom 
Ij*  0 


S  (5fV)  ♦  ^P-(2p^) 

u  • 

E.  •  6.11  eV 
(63  acates) 


‘Lu(5£37a6d)  ^Pj<2p^) 

E.  -  6.146  eV 
(81  atates) 


S^(5£57a6d)  +  ^j(2p*) 

b  •  6.223  eV 
(72  atates) 


Molecular  States 


V(l).  V(2),  V(l),  V(2).  V(l). 

V(2),  ^n(3),  *n(3),  *n(3),  ^a(3), 
*6(3).  ^6(3),  *6(3) .  S(3),  ^r(3),  *r(3) 
‘r(3).  ^(3),  *H(3).  ‘h(3),  ^1(2),  S(2) 
*1(2),  ^(l).  *K(1),  ^K(l) 

^r(l),  *1*^(2),  S'(l),  *r'*“(2),  V(l), 

*1I'3).  *11(3).  ®n(3),  *A(3),  *A(3),  ®d(3) 
*♦(3).  S3),  S3).  *r(3),  *r(3).  ®r(3) 
*H(3).  *H(3).  ®H(3).  *1(3),  *1(3),  ®1(3) 
*K(3),  *K(3).  ®K(3),  *A(2),  *A(2),  ®A(2) 
Sl(l),  *M(1).  ®M(1) 

*r-(2),  V(i).  *r’(2),  *1*^(1) .  ®r"(2), 

*n(3).  *0(3).  ®n(3).  *£(3).  *A(3),  ®6(3) 
S(3),  *6{3),  ®6(3).  *r(3>,  *r(3),  ®r(3) 

*H(3),  ®i(3).  ®H(3),  *1(3),  *1(3),  ®I(3) 
*K(2),  *K(2),  *fc(2),  *A(1),  *A(1),  ®A(1) 


Table  2*  Molecular  states  of  UO"*"  (Concluded) 


Separated  Atom  Molecular  States 

ir^  +  o~ 

5Lu(5f36d)  +  2p^(2p5)  V(l),  V(2),  ^r+d),  h- (2) ,  ^11 

(54  states)  ^n(3),  ^A(3),  ^A(3),  ^*(3),  ^«(3), 

^r(3),  hO),  ^H(3),  6h(3),  ^1(3), 
6i(3),  ^K(3),  6k(3).  ^A(2),  ^h(2) , 

S<(1),  6m(1) 

+  0” 


^Iu(5f3)  +  ^g(2p^)  H-ii),  ^n(i),  ^A(i),  ^*(1).  ^r(i) 

(8  states)  ^H(l),  ^I(l),  ^K(l) 


(1^  ( cm”  ^ ) : 

a)gXg(cin“^): 

R,(A): 


Table  3.  Calculated  spectropic  data  for  UO/UO'*’ 


820  (matrix) 

836  (scaled  from  ThO) 
863  (SCF)* 

845  (MCSCF)** 

859  (RDF-UTRC) 

2.7  (H-H  potential) 
12.0  ( excited  state) 


UQ-*-  [‘*1 

949  (SCF) 

935  (MCSCF) 

925  (MCSCF:  Q  =  9/2) 
980  (RDF-UTRC) 


2.7  (H-H  potential) 


1.84  (matrix)  1.83  (SCF) 

1.838  (Heaven  and  Nicolai)** 

1.88  (SCF)  1.84  (MCSCF) 

1.84  (MCSCF-Krauss)  1.84  (RDF-UTRC) 

1.89  (RDF-UTRC) 


*SCF  and  MCSCF  results,  see  Reference  78. 


Table  4.  5f  spin-orbit  parameters  for  U  ions 


Conf igurat ion 

4  (cm 

LASL  (Wood,  Boring) 

(5f) 

2442 

(5f2) 

2212 

(5f3) 

1977 

C  =  — ^  ); 

;  k*  =  A  for  j  -  1/2 

2U1 

k  =  -(i-  +  1)  for  j  +  1/2 

Reference  53 


Table  5.  Terra  levels  for 

U+3  (5f  ),  cm 

Terra 

Energy 

2f 

105180. 

2g 

72366, 

2f 

58586, 

2d 

48896, 

% 

47659. 

2l 

42749. 

“♦D 

42566. 

2l 

41939. 

2p 

32066. 

2d 

32060. 

2k 

26419. 

“♦G 

26097. 

2g 

23665. 

2h 

18222. 

16782. 

“♦F 

16653. 

0. 

=  33343  cra-1. 

avg 
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Table  6.  Excitation 


Energies  of  U’*’**» 


Configuration 


E  (cm 


1983  cm“M 
2873  cm"^) 


5£S/2 


5f^5/2  ^^7/2 
5^5/2  5^^7/2 
5f^/2 

5f5/2^^7/2^‘*3/2 
5f5/25f7/2^‘*5/2 
5f 27/2683/2 
5f27/26d5/2 


6938. 

13898. 

28890. 

30592. 

znik. 

38166. 

45342. 

45754. 

52938. 


(  =  2220  cm~^) 

<=60  -  ”53  ■="-') 


5' ‘5/2 

5^5/2^^7/2 

5£S/2 

515/2683/2 

5f3/26d5/2 

517/2683/2 

517/2685/2 


15570. 

60032. 

68410. 

68912. 

77302. 


2455  cm” 
4268  cm" 


^^5/2 

5f7/2 

683/2 


92921. 

103590. 
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Table  7.  Term  levels  for  central  atom  (one-electron)  excitation  of  UO^ 


Molecular  Orbital  Symmetry  Te  (cm~^) 


1  <))U 

0. 

1  6u 

-362.* 

3mi 

'"u 

10118. 

4og 

32250. 

16g 

33163. 

1  qg 

40976. 

2TTg 

'"g 

45695. 

*  Although  the  center-of-gravity  of  the  1 6u  state  lies  lower  than  l(l)u,  spin- 

2 

orbit  splitting  causes  the  become  the  ground  state,  in 

agreement  with  previous  studies  of  this  system  (Reference  81) 


fc  ^ 

k 


> 

% 

$ 

•\ 

■.'* 

K' 

l; 


Table  9.  Calculated  assignments  of  electronic  transitions  in  UO2* 


Transition 


"♦5/2u  " 

^3/2u 

1 

^*5/2u  " 

^^5/2u 

1' 

^*5/2u  ■*■ 

^*7/2u 

i 

I 

2*  - 
*5/2u  ^ 

^"3/2u 

2*  . 
*5/2u  ^ 

^^3/2g 

1 

2*  + 

*5/2u 

^"3/2g 

(f  +  d) 


<‘-g  Ha  ♦  t'l  hJ) 

2»  ^  Av| 

*5/2  "■  «5/2g 
'*'g  ‘*0  *  ‘i  H„3'u) 


Wavelength 


Strength 


>  5u  (calc) 

1.5  -  1.611  (calc) 
l,45ii  (calc) 

■“1.511  (exp) 

0.84p  (calc) 

■“0.911  (exp) 

0.3011  (calc)  Moderately  Strong  (u  •*•  g) 

Strong  (u  g) 


Weak  (u  +  u) 
Weak  (u  •♦•  u) 
Weak  (u  ♦  u) 

Weak  (u  ■♦■  u) 


0.18p  (calc) 


O.I611  (calc) 


Very  Strong  (u  +  g) 
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Calculated  spectroscopic  data  for  UO 

ojj  (cni~^)  885.4 

Xjj  (cm~^)  -2.7 

Fg  (A)  1.73 

D  (cm“^)  62100.  (exp) 
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Table  11.  Calculated  oscillator  strengths  for  the  vibrational  transitions  of  U02-^  (Concluded) 
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Table  12.  Total  integrated  absorption  coefficents  for  the  UOj'*’  ground  state 


Absorption  Coefficient,  cm" 2  atm” ^ 


Temperature , 
(Wavelength, 

*K 

JL 

Fundamental 

(11.43) 

First 

Overtone 

(5.73) 

Second 

Overtone 

(3.83) 

Third 

Overtone 

(2.88) 

Total 

100. 

2773.12 

8.51 

0.05 

0.00 

2781.68 

273.15 

2112.97 

8.68 

0.05 

0.00 

2781.70 

300. 

2772.89 

8.77 

0.06 

0.00 

2781.71 

500. 

2771.71 

10.01 

0.07 

0.00 

2781.79 

1000. 

2766.50 

15.31 

0.14 

0.00 

2781.95 

1500. 

2760.16 

21.61 

0.27 

0.00 

2782.04 

2000. 

2753.40 

28.23 

0.45 

0.01 

2782.09 

2500. 

2746.39 

35.01 

0.69 

0.02 

2782.10 

3000. 

2739.19 

41.90 

0.99 

0.03 

2745.12 

4000. 

2724.32 

55.95 

1.78 

0.07 

2782.12 

5000. 

2708.17 

70.16 

2.81 

0.15 

2781.29 
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Field  Command,  Defense  Nuclear  Agency 

ATTN 

Code  F31 

ATTN 

FCPR 

ATTN 

Code  X211,  Tech  Library 

ATTN 

FCTT 

ATTN 

FCTT,  W.  Surma 

Space 

S  Naval  Warfare  Systems  Cmd 

ATTN 

PME  117-20 

Joint  Strat  Tgt  Planning  Staff 

ATTN 

JLKS 

DEPARTMENT  OF  THE  AIR  FORCE 

Under  Secy  of  Def  for  Rsch  &  Engrg 
ATTN:  Defensive  Systems 
ATTN:  Strat  &  Space  Sys  (OS) 

ATTN:  Strat  &  Theater  Nuc  For,  F.  Vajda 

DEPARTMENT  OF  THE  ARMY 


Atmospheric  Sciences  Laboratory 
3  cys  ATTN:  OELAS-EO,  F.  Niles 

BMD  Advanced  Technology  Center 
ATTN:  ATC-O,  W.  Davies 
ATTN:  BMDATC-D,  M.  Capps 


Force  Geopivys 

ATTU.  I  CD 


ATTN 
:ys  ATTN 
;ys  ATTN 
:ys  ATTN 
:ys  ATTN 
:ys  ATTN 


LSP 

LID 

LKO 

LS. 

LSI, 

LSP, 


lies 

D. 


Laboratory 

Smith 

Swider 

Huffman 

Murphy 

Sharma 

Nadile 

Champion 

Paulson 

Stair 


Air  Force  Office  of  Scientific  Rsch 
ATTN:  AFOSR/NC 
ATTN:  AFOSRINP 


BMD  Systems  Cotmiand 
2  cys  ATTN:  BMDSC-HW 

Dep  Ch  of  Staff  for  Rsch  Dev  &  Acq 
ATTN:  DAMA-CSS-N 

Harry  Diamond  Laboratories 
2  cys  ATTN:  DELHD-NW-P 

US  Army  Ballistic  Research  Lab 

ATTN:  DRDAR-BLA-S,  Tech  Library 
ATTN:  DRDAR-BLB,  J.  Mester 

US  Army  Combat  Surv  &  Target  Acq  Lab 
ATTN:  DELCS-K,  S.  Kronenberg 

US  Army  Foreign  Science  &  Tech  Ctr 
ATTN:  DRXST-SD 

US  Army  Nuc  ear  &  Chemical  Agency 
ATTN:  Library 


Air  Force  Systems  Command 
ATTN :  DLAE 
ATTN :  DLS 
ATTN:  DLTW 
ATTN:  DLXP 
ATTN:  SDR 

Air  Force  Technical  Applications  Ctr 
ATTN:  TD 

ATTN:  Tech  Library 
ATTN:  TF 

Air  Force  Weapons  Laboratory 
ATTN:  SUL 

Air  University  Library 
ATTN:  AUL-LSE 

Dep  Ch  of  Staff,  Rsch,  Dev,  &  Acq 
3  cys  ATTN:  AFRDS,  Space  Sys  &  C3  Dir 


US  Army  Research  Office 
ATTN:  R.  Mace 


Rome  Air  Development  Center 
ATTN:  OCD,  J.  Simons 
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Space  Division/AQ 
ATTN;  WE 

Space  Division 

ATTN:  YGD 
ATTN;  YN 

Strategic  Air  Conmand 
ATTN:  SAC/S 10 
ATTN:  INCR 
ATTN:  NRI/STINFO 
ATTN:  XPFR 
ATTN:  XPFS 

DEPARTMENT  OF  ENERGY 


OTHER  GOVERNMENT  AGENCIES  (Continued) 

Dep  of  Conmerce,  Natl  Oceanic  &  Atmospheric  Admin 
3  cys  ATTN:  E.  Ferguson 
3  cys  ATTN:  F.  Fehsenfeld 

Institute  for  Telecommunication^  Sciences 
ATTN:  G.  Falcon 
ATTN:  W.  Utlaut 

NASA,  George  C.  Marshall  Space  Flight  Center 
ATTN:  0.  Watts 
ATTN:  N.  Stone 
ATTN:  W.  Roberts 

NASA,  Ames  Research  Center 

ATTN:  N-245-3,  R.  Whitten 


Department  of  Energy 

Office  of  Military  Application,  GTN 
ATTN:  OMA,  DP-22 

University  of  California 

Lawrence  Livermore  National  Laboratory 
ATTN;  L-10,  A.  Grossman 
ATTN:  L-10,  H.  Kruger 
ATTN;  L-262,  D.  Wuebbles 
ATTN:  L-325,  G.  Haugan 
ATTN:  L-48,  E.  Woodward 
ATTN:  L-71,  J.  Chang 

Los  Alamos  National  Laboratory 
ATTN:  D.  Sappenfield 
ATTN:  G.  Smith 
ATTN:  Librarian 
ATTN:  M.  Pongratz 
ATTN:  M.  Sanford 
ATTN:  MS362  Library 
ATTN:  R.  Jeffries 
ATTN;  T.  Bieniewski 
ATTN:  T.  Kunkle,  ESS-5 

Sandia  National  Laboratories 
ATTN:  T.  Cook 

Sandia  National  Laboratories 
ATTN:  L.  Anderson 
ATTN:  M.  Kraimi 
ATTN:  Org  1250,  W.  Brown 
ATI,  Tech  Library  3141 

OTHER  GOVERNMENT  AGENCIES 

Central  Intell igence  Agency 
ATTN:  OSWR/NED 


NASA,  GSFC/WFF 

ATTN:  J.  Gray 

NASA  Headquarters 

ATTN:  Schardt,  Code  EE 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aero-Chem  Research  Labs,  Inc 
ATTN:  A.  Fontijn 

Aerodyne  Research,  Inc 
ATTN:  C.  Kolb 
ATTN:  J.  Wormhoudt 
ATTN:  M.  Camac 
ATTN:  M.  Zahnhitzer 

Aerospace  Corp 

ATTN:  H.  Mayer 
ATTN:  I.  Garfunkel 
ATTN:  J.  Reinheimer 
ATTN:  J.  Straus 
ATTN:  V.  Josephson 
ATTN:  N.  Cohen 

AVCO  Everett  Research  Lab,  Inc 
ATTN:  A830 

ATTN;  C.  von  Rosenberg,  Jr 
ATTN:  Tech  Library 

Battel le  Memorial  Institute 
ATTN:  R.  Thatcher 
ATTN:  STOIAC 

Berkeley  Rsch  Associates,  Inc 
ATTN:  J.  Workman 
ATTN:  S.  Brecht 


Dep  of  Commerce,  National  Bureau  of  Standards 
ATTN:  A.  Phelps 

Dep  of  Conmerce,  National  Bureau  of  Standards 


ATTN: 

Sec 

Ofc 

for  J. 

Devoe 

ATTN: 

Sec 

Ofc 

for  M. 

Krauss 

ATTN: 

Sec 

('  ‘■c 

for  R. 

Levine 

ATTN: 

Sec 

Ofc 

for  S. 

Abramowi tz 

NASA,  Goddard  Space  Flight  Center 


ATTN 

ATTN 

ATTN 


Code  6801,  A.  Tempkin 
Code  900,  J.  Siry 
Tech  Library 
a  -  • 


Boston  College,  Space  Data  Analysis  Lab 
ATTN:  E.  Hegblom 

ATTN:  W.  Grieder 

The  Trustees  of  Boston  College 
2  cys  ATTN:  Chairman  Dept  of  Chemistry 

2  cys  ATTN:  Chairman  Dept  of  Physics 

University  of  California  at  Riverside 
ATTN:  J.  Pitts,  Jr 


California  Institute  of  Technology 
ATTN:  J.  Ajello 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


Calspan  Corp 
ATTN: 
ATTN: 
ATTN: 
ATTN: 


C.  Treanor 
J.  Grace 
M.  Dunn 
W.  Wurster 


Mission  Research  Corp 


University  of  Colorado 

ATTN:  C.  Lineberger,  JILA 
ATTN:  G.  Lawrence,  LASP 

Computer  Sciences  Corp 

ATTN:  F.  Eisenbarth 

Concord  Sciences  Corp 
ATTN:  E.  Sutton 

Cornell  University 

ATTN:  M.  Kelly 

University  of  Denver 

ATTN:  Sec  Officer  for  D.  Murcray 

University  of  Denver 

ATTN:  B.  van  Zyl 

Environmental  Rsch  Inst  of  Michigan 
ATTN:  IRIA  Library 

EOS  Technologies,  Inc 

ATTN:  B.  Gabbard 

General  Electric  Co 

ATTN:  P.  Zavitsanos 
ATTN:  R.  Edsall 

Geo  Centers,  Inc 

ATTN:  E.  Marram 

HSS,  Inc 

ATTN:  H.  Shuler 

Institute  for  Defense  Analyses 
ATTN:  E.  Bauer 
ATTN:  H.  Wolfhard 

IRT  Corp 

ATTN:  H.  Mitchell 


ATTN:  H.  Dickinson 

Kaman  Sciences  Corp 

ATTN:  E.  Conrad 

Kaman  Tempo 

ATTN:  B.  Gambill 
5  cys  ATTN:  DAS  I AC 

Kaman  Tempo 

ATTN:  DAS  I AC 

Lockheed  Missiles  &  Space  Co,  Inc 
ATTN:  B.  McCormac 
ATTN:  J.  Cladis 
ATTN:  J.  Kumer 
ATTN:  J.  Perez 
ATTN:  J.  Reagan 
ATTN:  M.  Walt 
ATTN:  R.  Sears 


ATTN 
ATTN 
ATTN 
ATTN 
ATTN 
ATTN 
ATTN 
2  cys  ATTN 


D.  Archer 
D.  Sowle 
F.  Guigliano 
M.  Scheibe 
P.  Fischer 
R.  Hendrick 
R.  Kilb 
Tech  Library 


Mission  Research  Corp 


C.  Longmire 
R.  Peterkin 
R.  Stellingwerf 


Nichols  Research  Corp,  Inc 
ATTN :  R .  Byrns 

Pacific-Sierra  Research  Corp 

ATTN:  H.  Brode,  Chaiman  SAGE 

Photometries,  Inc 

ATTN:  I.  Kofsky 

Physical  Dynamics,  Inc 

ATTN:  A.  Thompson 

Physical  Research,  Inc 

ATTN:  T.  Stephens 

Physical  Research,  Inc 
ATTN:  J.  Devore 

Physical  Science  Lab 

ATTN:  W.  Berning 

Physical  Sciences,  Inc 
ATTN:  D.  Greem 
ATTN:  G.  Caledonia 
ATTN:  K.  Wray 
ATTN:  R.  Taylor 
ATTN:  T.  Rawlings 

University  of  the  Coinnonwealth,  Pittsburgh 
ATTN:  F.  Kaufman 
ATTN:  M.  Biondi 
ATTN:  W.  Fite 

Princeton  University 

ATTN:  Librarian 


RSD  Associates 


F.  Gilmore 
H.  Dry 
R.  Lindgren 
R.  Turco 


R&D  Associates 

ATTN:  B.  Yoon 
ATTN:  J.  Rosengren 

Rand  Corp 

ATTN:  C.  Crain 
ATTN:  P.  Davis 

Rand  Corp 

ATTN:  B.  Bennett 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Science  Applications  Inti  Corp 
ATTN:  D.  Kara! in 

Science  Applications  Inti  Corp 
ATTN:  E.  Hyman 

Space  Data  Corp 

ATTN:  S.  Fisher 

SRI  International 

ATTN:  C.  Rino 
ATTN:  D.  McDaniel 
ATTN:  J.  Casper 
ATTN:  J.  Vickrey 
ATTN:  President 
ATTN:  R.  Leadabrand 
ATTN:  W.  Chesnut 

SRI  International 

ATTN:  C.  Hulburt 

Stewart  Radiance  Laboratory 
ATTN:  R.  Huppi 

Technology  International  Corp 
ATTN:  W.  Boquist 

Teledyne  Brown  Engineering 
ATTN:  F.  Leopard 
ATTN:  MS-12  Tech  Library 
ATTN:  N.  Passino 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Toyon  Research  Corp 
ATTN:  J.  Ise 

United  Technologies  Research  Ctr 
2  cys  ATTN:  H.  Michels 
2  cys  ATTN:  R.  Hobbs 


Utah  State  University 


ATTN 

A.  Steed 

ATTN 

C.  Wyatt 

ATTN 

D.  Baker 

ATTN 

K.  Baker,  Dir  Atmos  &  Space  Sci 

Vi  Si  Dyne,  Inc 

ATTN 

J.  Carpenter 

Wayne  State  University 

ATTN 

R.  Kuiranler 

Wayne  State 

University 

ATTN 

W.  Kauppila 

DIRECTORY  OF  OTHERS 


Government  Publications  Library-M 
ATTN:  J.  Winkler 

Yale  University 

ATTN:  Engineering  Department 
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